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Abstract: Within the gauge/gravity correspondence, we discuss the general formulation of 
the shockvi^ave metric which is dual to a 'nucleus' described by the strongly-coupled = 4 
SYM theory in the limit where the number of colors Nc is arbitrarily large. We emphasize that 
the 'nucleus' must possess degrees of freedom per unit volume, so like a finite-temperature 
plasma, in order for a supergravity description to exist. We critically reassess previous proposals 
for introducing transverse inhomogeneity in the shockwave and formulate a new proposal in 
that sense, which involves no external source but requires the introduction of an 'infrared' cutoff 
which mimics confinement. This cutoff however plays no role when the shockwave is probed by a 
highly virtual projectile, so like in deep inelastic scattering. We consider two such projectiles, the 
dilaton and the 7^-current, and compute the respective structure functions including unitarity 
corrections. We find that there are no leading-twist contributions to the structure functions 
at high virtuality, meaning that there are no point-like constituents in the strongly coupled 
'nucleus'. In the black-disk regime at low virtuality, the structure functions are suggestive of 
parton saturation with occupation numbers of order one. The saturation momentum Qs grows 
with the energy like Q'^ 1/x (with x the Bjorken variable), which is the hallmark of graviton 
exchanges and is also necessary for the fulfillment of the energy-momentum sum rules. 
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1. Introduction 

Some of the experimental discoveries at RHIC, notably the unexpectedly large medium effects 
known as elliptic flow and jet quenching, led to the suggestion that the deconfined hadronic 
matter produced in the intermediate stages of a heavy ion collision might be strongly interact- 
ing [1,2]. This observation triggered a large theoretical activity aiming at understanding the 
properties of a strongly-coupled plasma via string theory methods, within the AdS/CFT corre- 
spondence (see, e.g., the review papers [3-5] and refs. therein). Part of this activity concentrated 
on the 'dual' formulation of a heavy ion collision at strong coupling. The main motivation for 
that problem is less in the description of the nuclear wavefunctions by themselves — a real 
ultrarelativistic nucleus in QCD is rather described via weak coupling methods, because of the 
high parton density produced by the high-energy evolution [6] — , but rather in understanding 
the late stages of the collision and, especially, study thermalization. If the deconfined matter 
produced after a heavy ion collision is effectively strongly coupled, then one may hope that the 
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subsequent dynamics leading to thermalization is not very sensitive to the perturbative early 
stages, and thus it can be studied by simulating a full collision in AdS/CFT. Such a study may 
help explaining one of the main puzzles left by RHIC, which is the very short thermalization 
time suggested by the elliptic flow data. 

The AdS/CFT correspondence [7-9] deals with the conformally symmetric = 4 super- 
symmetric Yang-Mills (SYM) theory, and not directly with QCD. But thus is probably not a 
major impediment so long as we are interested in the temperature range 2Tc ^ T < 5Tc (the 
relevant range for heavy ion experiments at RHIC and LHC), where the QCD plasma itself 
appears to be nearly conformal, as shown by lattice studies [10]. For practical reasons, all the 
studies so far were limited to the 'supergravity approximation' valid when the number of colors 
Nc, and hence also the 't Hooft coupling A = g^Nc, are arbitrarily large. (The gauge coupling 
g is assumed to be small, g 1^ but fixed.) In this approximation, a plasma or a nucleus are 
described as modifications of the asymptotically AdS5 geometry that the string theory lives in, 
and the dynamics reduces to classical gravity in this curved space-time. In order for the re- 
spective metric perturbation to survive as A'c ^ oo, one needs to assume that the system under 
consideration possesses A^^ degrees of freedom per unit volume. This is indeed the case for the 
= 4 SYM plasma at finite temperature, whose dual description is a black-brane in AdSs [11], 
but it would not be true for a QCD-like nucleus, which is built with colorless 'hadrons'. 

Still, if we are merely interested in studying thermalization, then one can replace the colliding 
nuclei with finite-size slices of plasma, represented as slices of 'black-holes' (more properly, 
black branes) in AdSs. In the center ~of-mass frame, where these 'nuclei' (prior to the collision) 
are highly energetic and hence strongly Lorentz contracted, the black-hole slices appear as 
Shockwaves directing against each other. The ultimate objective is to explicitly compute the 
scattering between two such Shockwaves by solving Einstein equations, and thus investigate the 
evolution towards a black brane at late stages, which is synonymous of thermalization. Several 
steps have been already done in that sense, with interesting results, and the dynamical formation 
of a black hole has been already seen within supergravity calculations [12-21]. 

But before attacking this complicated problem, it is the very formulation of the gravity dual 
of a nucleus as a shockwave which still meets with difficulties in the literature. The would-be 
natural recipe for building such a shockwave in view of the discussion above — namely, start 
with a slice of the plasma with finite longitudinal width L ^ 1/T in the plasma rest frame and 
then boost the dual 'black-hole slice' metric up to a large Lorentz factor 7 S> 1 — is not fully 
satisfactory at a conceptual level, since a 'slice of a black-hole' in AdSs is not an exact solution 
to Einstein equations. Yet, as we shall demonstrate in Appendix this approximate procedure 
has the merit to generate a shockwave metric which is similar to the one that we shall propose 
via different considerations, in Sect. 2. 

An important related problem is that of the conformal symmetry breaking: this was already 
necessary in order to built hadronic-like bound states in the = 4 SYM theory [22-26], and 
one can hardly see how one could construct a nuclear shockwave without introducing any scale 
in the problem. Yet, the original proposal in that sense, due to Janik and Peschanski [12], seems 
not to require any such a scale. As we shall argue in Sect. 2, such a scale must nevertheless be 
added by hand to that construction, in the form of a cutoff on the radial dimension^ of AdS5 — 

^Throughout this paper, we will be using Fefferman-Graham coordinates (see Sect. 2 for precise definitions), 
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either a 'hard wah', i.e. a sharp radial cutoff at z = 1/A [23], or a 'soft wall', like in Rcf. [25]. 
This cutoff corresponds to an infrared cutoff A in the M = 4 SYM theory and mimics the 
confinement physics of QCD. Incidentally, the metric produced by boosting a slice of the plasma 
is identical to that by Janik and Peschanski near the Minkowski boundary (for z ^ l/T), but it 
naturally involves a radial cutoff, namely the black-hole horizon at z ^ 1/T (see Appendix A). 

The need for a radial cutoff in AdSs can be also seen via a different argument, which is 
internal to the supergravity construction: without such a cutoff, the shockwave metric develops 
a genuine space-time singularity at z ^ oo, and not just a coordinate singularity. This has been 
recognized in the literature [27] for the analog of the Janik-Peschanski metric in 4 dimensions, 
known as the Kaigorodov space-time [28], but a similar argument holds in 5 dimensions. The 
argument is quite subtle, since the singularity at 2; — > oo does not show up in any of the curvature 
scalars computed from this metric, but is of a more general type, known as ^-incompleteness [29]. 

The issue of the radial cutoff becomes even sharper when trying to introduce transverse 
inhomogeneity in the shockwave, corresponding to a profile for the nuclear matter distribution 
in the plane transverse to the collision axis. The two types of proposals in that sense that we 
are aware of [14-17, 30] use external sources to break down conformal symmetry. Besides the 
difficulty to motivate the iV^-scaling for the strength of these sources, these proposals have some 
unwanted features, which make them unrealistic from the viewpoint of QCD. The proposal in 
Ref. [14] uses a pointlike source located in the bulk of AdS5 to produce a 'nuclear' energy- 
momentum distribution which falls off at large transverse distances according to a power law, 
instead of the exponential law expected in the presence of confinement. Accordingly, when 
used for scattering problems, this shockwave yields cross-sections which violate the Proissart 
bound. The proposal in Ref. [30] is even more objectionable, in that one is mistreating the 
'ultraviolet' (high-momentum) part of the 'nuclear' wavefunction. This proposal involves a 
source on the boundary of AdSs, hence the corresponding metric perturbation docs not vanish 
when approaching the boundary. Via the UV/IR correspondence [26,31-35], this means that 
the spectrum of the quantum modes included in the dual 'nuclear' wavefunction is flat in the 
high-momentum limit, rather than rapidly falling down, as it should on physical grounds. 

On the other hand, the analysis in Ref. [30] has the virtue to have identified a class of exact 
solutions to Einstein equations of the shockwave type, which allow for a generic inhomogeneity in 
the transverse plane. This analysis lies at the basis of our shockwave proposal in this paper. One 
can succinctly describe our proposal as follows: among the two general solutions found in [30] 
and which involves the modified Bessel functions K2 and, respectively, I2, we shall discard the 
solution oc K2 that was adopted in Ref. [30] and keep the other solution oc I2. The latter is the 
generalization of the shockwave metric by Janik and Peschanski [12] (which is homogeneous in 
the transverse plane) to a generic transverse energy-momentum distribution. This solution van- 
ishes near the boundary (z 0) like z^, meaning that the spectrum of the modes included in the 
dual 'nucleus' has an acceptable high-momentum tail oc On the other hand, this solution 

would exponentially blow up at large values of z, but this is not a problem since, as previously 
argued, the radial dimension must be anyway supplemented with a cutoff at 2; = 1/A. From the 
perspective of the boundary gauge theory, this scale A plays several roles, so like in QCD: it acts 

with the radial distance denoted as z. Thus by 'large radial distances' we mean large separations from the 
Minkowski boundary of AdSs, which in these coordinates lies at 2: = 0. 
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as an infrared cutoff for quantum fluctuations, it sets the scale for energy-momentum density 
in the 'nucleus', and also the characteristic scale for transverse inhomogcneity. 

The introduction of this radial cutoff is clearly an ad hoc procedure (it spoils the exactness 
of the solution) and any quantity which is sensitive to the details of this procedure is model- 
dependent. This is probably the case for the collision between two such Shockwaves, but this is 
hardly a surprise: already in QCD, total hadronic cross-sections are dominated by soft interac- 
tions and thus are sensitive to the physics of confinement. But even in that case, one may hope 
that the thermalization process in the late stages of the collision is less sensitive to the details 
of the cutoff. On the other hand, the collision between this 'nucleus' and a 'hard' probe with 
high transverse resolution Q S> A is probing the Shockwave metric only at small radial distances 
-^^^l/Q^l/-^) S'lid hence it is completely insensitive to the infrared cutoff. Once again, this 
is similar to QCD where deep inelastic scattering at high-Q^ is insensitive to confinement. 

The above considerations arc illustrated by our calculations of deep inelastic scattering 
(DIS) off the Shockwave, which represent most of the material in this paper. The first calcula- 
tions of DIS at strong coupling within the gauge/string duality referred to some other types of 
targets: a 'glueball' bound state [33, 36-38] (whose gravity dual is a normalizable dilaton state 
in AdSs) and a finite temperature plasma (as represented by the AdSs black-brane) [34,35]. 
More recently, Refs. [39,40] presented the corresponding calculations for a shockwavc target (a 
'nucleus', or a 'plasma slice'), but their respective results appear to disagree with each other. 
The main difference refers to the energy dependence of the saturation momentum Q.,, which 
is the characteristic scale for the onset of unitarity corrections: for Q ^ Qs, the scattering is 
weak, whereas for Q < Qs it reaches the unitarity bound, or 'black disk', limit. 

Already at weak coupling, Qg is known to grow quite fast with the energy, as an inverse 
power of the Bjorken x variable, due to the rapid increase in the gluon density at small x via 
bremsstrahlung [6]. At strong coupling, one expects parton branching to be faster and to 'quasi- 
democratically' divide the energy among the daughter partons [33,34,36,41-43]. Accordingly, 
all partons should fragment down to very small values of x and then the saturation momentum 
should grow with 1/x even faster. This is indeed what one found in Ref. [34,36,40] for various 
types of hadronic targets: dilaton, infinite plasma, and a plasma slice. Moreover, the energy 
dependence emerging from these calculations looks very natural from the viewpoint of super- 
gravity: the scattering proceeds via multiple graviton exchanges, so the saturation momentum 
inherits the energy dependence of the graviton propagator, which implies^ oc 1/x. This 
specific energy dependence is also necessary to ensure energy-momentum conservation, as we 
shall later explain. On the other hand, Ref. [39] reported a rather counterintuitive result for Qg, 
which becomes independent of x in the high energy limit. Although that analysis uses a dif- 
ferent projectile — namely, a small 'color dipole' (dual to a Nambu-Goto open string in AdSs) 
instead of a virtual photon — , it is unlikely that this can explain the dramatically different 
result obtained there for the saturation momentum^. Indeed, the latter is an intrinsic property 

■^The picture is more subtle in the case of an infinite plasma, and the associated saturation scale shows a faster 
rise with the energy, namely Ql oc l/x^, but the additional power of 1/x is understood simply as the coherence 
length of the virtual photon [34, 35] . 

■'in fact, we shall find that the supergravity field dual to the virtual photon describes a dipolar partonic 
fluctuation of the latter, so the DIS proceeds via the dipole scattering off the Shockwave, so like in QCD at weak 
coupling. Hence even the physical setup looks a priori similar in Refs. [34,36,40] and respectively [39]. 



-4- 



of the target, hence it must be insensitive to the precise nature of the projectile. 

Our analysis in Sects. 3 and 4 will extend the previous analysis [40] of a finite-width plasma, 
with results which agree with Ref. [40] whenever a direct comparison is possible. Our general- 
ization will refer to several aspects: 

(i) We shall consider two types of external probes: an 7?.-current (so like in Ref. [40]) in 
Sect. 4, and a 'dilaton' in Sect. 3. The 7^~current is an analog of the electromagnetic current for 
the = 4 SYM theory, and is dual to an Abelian vector field propagating in the asymptotically 
AdSs space-time. In the supergravity approximation, DIS amounts to solving the Maxwell 
equations for this vector field in the background of the AdSs shockwave geometry'*. The 'dilaton' 
is a supergravity scalar field which is dual to the Lagrangian density C = (Xj^Flj^^Fa^ + . . . . 
The dilaton case will be presented first, and in more detail, since the corresponding equations 
of motion — the Klein-Gordon equations in the shockwave geometry — turn out to be simpler. 

(ii) We shall for the first time consider the impact parameter dependence of the scattering 
amplitude and of the saturation momentum, and in particular study the expansion of the black 
disk with increasing energy and the emergence of the Froissart bound. We shall perform this 
analysis for two types of Shockwaves: that introduced in Ref. [14] and which involves a pointlike 
source in the bulk, and the one that we shall propose in Sect. 2 below and which allows for an 
arbitrary inhomogeneity in the transverse plane. 

(iii) Both for the dilaton and for the 7^-current, we shall construct the exact supergravity 
solution in great detail, by resumming multiple scattering to all orders and showing that this 
leads to the eikonalization of the single graviton exchange, as expected for the scattering off 
a shockwave. (The eikonalization was also advocated in Ref. [40], but without an explicit 
construction.) This requires the respective 'bulk-to-bulk' propagators in AdSs in real time and 
light -cone coordinates, that will be constructed in Appendix^. Similar results in the AdS/CFT 
context have been previously obtained in Refs. [44-46]. 

Let us now summarize the main results and conclusions which will emerge from our analysis: 
The DIS structure functions in the supergravity approximation come out in a factorized 
structure which is reminiscent of the A;_l -factorization in QCD at weak coupling [6], with the 
radial distance z in AdSs playing the role of the transverse size r_|_ ~ of the partonic 

fluctuations of the projectile ('virtual photon') in the boundary gauge theory (as expected from 
the UV/IR correspondence [26,34,35]). Namely, they involve a convolution in z over the square 
of the incoming field (the supergravity solution for the dilaton or the Maxwell field in the absence 
of the shockwave) times the cross-section for the scattering between this field and the shockwave 
at a given value oi z. In turn, this cross-section is expressed in the eikonal approximation, as 
the integral over all the impact parameters of a scattering amplitude which resums multiple 
graviton exchanges to all orders. 

For sufficiently large and/or large values of Bjorken x, the multiple scattering series can 
be expanded out, and then the dominant contribution to the structure functions comes from 
the double graviton exchange (the single graviton exchange being purely real). This leading 
contribution is of higher twist order, meaning that there are no point-like constituents in the 
target, in agreement with previous analyses at strong coupling [33-38], but in sharp contrast 

''For the supergravity approximation to apply to DIS, the total COM energy squared s must be limited in 
such a way not to allow for the excitation of massive string states. The respective condition can be written as 
^ -\/A [33], where x ~ Q^/s is the Bjorken variable. 
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with the situation at weak couphng, where the dominant contribution at large is of twist-two 
order and describes parton (in QCD, quark) distributions [6,47]. 

For sufficiently low at a given value of x, the amplitude reaches the unitarity limit and the 
structure functions are large. Remarkably, at least for the 7^-current, the structure functions at 
low have the same parametric form as for a proton in QCD at weak coupling. This similarity 
suggests a physical interpretation in terms of parton saturation at strong coupling [34,36] : the 
low-fcx and low-x region of the phase-space is filled with partons, with occupation numbers of 
order one. Note that, in pQCD, the occupation numbers at saturation are of order one only 
for quarks, but they are of order 1/A (with A = g^Nc <C 1 at weak coupling) for gluons [6]. In 
that case, saturation is driven by gluon dynamics, namely by the enhanced radiation of gluons 
with small values of x and their mutual interactions (quark saturate only due to their coupling 
to gluons). By contrast, at strong coupling, saturation is driven by quasi-democratic branching 
and the occupation numbers at saturation are of order one for all types of partons [34,36]. 

The borderline between the weak-scattering regime at high-Q^ and large-x, where there 
are no partons, and the saturation region at low-Q^ and low-x defines the saturation line, which 
can be expressed either as = Ql{x), or x = Xs{Q)- We find that this saturation line is exactly 
the same for the dilaton and the 7^-current, which confirms that this is an intrinsic property 
of the 'nuclear' target. Specifically, we find Ql{x) ~ A^L/x, in agreement with Ref. [40]; here, 
A is the 'confinement' scale, as introduced by the energy density in the target^, and L is the 
longitudinal extent of the target in its rest frame. This peculiar 1/x-dependence of Qs, which 
reflects the energy-dependence of the single graviton exchange, can be also understood via an 
independent argument, that we now explain. 

Namely, energy-momentum conservation implies that the integral dx F2{x , Q'^) of the 
structure function F2 has a flnite limit, of order A'^^, as — > cxd. At large x ^ Xg{Q), 
F2{x,Q'^) has only higher-twist contributions which rapidly die away with increasing (see 



Eq. ( 4.41 )). At X < Xs{Q), F2{x, Q^) is independent of x, due to saturation, and of order N^Q'^ 
(see Eq. ( [4.42D ). Hence, the integral is dominated by x ~ Xs{Q) — the highest value of x at 
which the structure function is still finite at large — , and can be estimated as 



/ dxF2(x,Q2) ^ XsF2{Xs,Q^) ~ XsN^Q"". (1.1) 

Jo 



For this to be independent of Q^, Xs must scale as Xs{Q) ~ or Q^(x) ~ 1/x, as announced. 

So far, we did not mention the dependence of the various results, so like Qs, upon the impact 
parameter b±. This will be discussed at length in the main text, and the main conclusion is 
that, at strong coupling, the function Q'^{x,b±) has the same 5_l -dependence as the energy- 
momentum distribution in the shockwave. In particular, for a distribution which exhibits an 
exponential tail at large b± (as expected in the presence of confinement), Q'^{x,b±) has an 
exponential tail as well, and the black disk area grows like ln'^(l/x), that is, it saturates the 
Froissart bound. On the other hand, for the shockwave metric in Ref. [14], Qs(x,6_l) ~ 1/6^ 
at large b±, and the area of the black disk grows like a power of the energy, in violation of the 
Froissart bound. 



*This is also the scale which fixes the radial cutoff at z ~ 1/A in AdSs, as previously explained, but this cutoff 
plays no role for DIS at > A^. 
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2. Shockwaves with impact parameter dependence in AdS^ 



In this section, we shall discuss various proposals for asymptotically AdS^ Shockwave metrics 
which are intended to represent the gravity duals of a fast moving 'nucleus' (more properly, a 

slice of the plasma) with a non-trivial profile in impact parameter space (the two dimensional 
space transverse to the direction of motion). We shall work in the nucleus infinite momentum 
frame, that is, we shall take the nucleus to move in the positive direction with a Lorentz 7 
factor which is arbitrarily large. It is then convenient to introduce light-cone coordinates, 

x+ = ^ (xO + x^), X- = ^ (x° - x^) , (2.1) 

in terms of which the nucleus moves in the positive x^ direction, and the only non-trivial 
component of its average energy-momentum tensor {T^^) is the (— ,— ) component (as this 
is the only one to be enhanced by the large factor 7^). By energy-momentum conservation, 

d^{Tfj,u) = 0, this component^ T = (T ) is independent of the light -cone time x"*". For the 

time being, we shall allow T [x^,x^) to be an arbitrary function of x^ and the transverse 

coordinates x_\_ = (x^,x^) (the 'impact parameter'). In practice we shall be mostly interested in 
the situation where the nucleus is strongly Lorentz contracted: T oc 5{x~) (the 'shockwave'). 

We shall assume T to be proportional to iV^. This is unrealistic from the point of view 

of a real QCD problem, where a nucleus is built with colorless hadrons (protons and nucleons). 
But within the present AdS/CFT context, such an assumption is necessary in order to be able to 
study the high-energy dynamics in the large- Ai'c limit, to which we would like restrict ourselves 

in what follows. Indeed, it is only when T oc A^^ that the metric perturbation induced by the 

'nucleus' in the 'bulk' metric of AdS5 is an effect of 0(1) which survives when Nc — 00. Then, 
the nucleus and its interactions can be simply described in the supergravity approximation, in 
which one first solves the Einstein equations to determine the asymptotically AdSs metric dual 
to the nucleus and then study the propagation of various projectiles in this metric. 

Our subsequent construction may be viewed as a model for the gravity dual of a slice 
of a = 4 SYM plasma at finite temperature T, which indeed has N"^ degrees of freedom 
per unit volume, and hence an energy density Too ~ A^^T^ in the plasma rest frame, yielding 

T ~ 7^A^r^ in the infinite momentum frame. In theoretical studies, the plasma is generally 

assumed to be infinite, but this is unrealistic from a phenomenological viewpoint and, besides, 
the concept of plasma makes sense also for a finite volume system (at least over a finite time 
interval), so long as the size L of the system along any direction obeys L ^ 1/T. An explicit 
connection between the shockwave metric that we shall construct in this section and the boosted 
slice of the plasma will be established in Appendix A. 

As mentioned in the Introduction, there are two types of proposals in the literature for 
shockwave metrics in AdSs : one which introduces a gravity source in the bulk of AdSs [14-17], 
and one which does not [12, 30] (but boundary sources are in principle allowed in the second 
case; see below). Both cases can be encoded in the following Einstein equations 

T^mn ~2 9mn '^^9mn ~ ^T^G^ Jjran i (^•^) 

^To avoid cumbersome notations, we shall omit the brackets denoting expectation values whenever there is no 
risk of confusion. 
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where m, n are 5~dimensional space-time indices, A = —Q/R^ (with R the curvature radius 
of AdSs) is the cosmological constant, TZ = g"^"'TZmn = — 20/i?^, G5 is the Newton constant 
in D = 5, and Jmn is the stress tensor of the source locaHzed in the bulk, for which we shall 
consider the two scenarios alluded to above: 

1. Type-I metric: Jmn = 0. This is the case considered in Rcfs. [12,30]. 

2. Type-II metric: corresponds to an ensemble of iV^ point-like 'particles' moving to- 
gether within AdSs along a null geodesic parallel to the Minkowski boundary (see Eq. ( 2.1^ ) 
for an explicit expression). This is the case considered in Refs. [14-17]. 

As shown in the literature, the corresponding shockwave solutions to Eq. ( |2.2| ) can be obtained 
with the following Ansatz (in the so-called Fefferman-Graham coordinates) 

ds^ = -^[dz'^ -2dx^dx- +dx]_ + h{z,x-,x±){dx'f], (2.3) 

for both cases: with or without bulk sources. Still in both cases, the function h{z,x~ ,x±) is 
determined by the following, linear, equation 

9^ - + Vi) /i(z,x-,xx) = -levrGsJ^-, (2.4) 



which is the (— ,— ) component of Eq. ( |2.2| ). Note that there is no 'small perturbation' as- 



sumption involved here: Eq. (2.4) is the exact consequence of the Einstein equations ( p.2| ) for 
the Ansatz ( ^^ ) and the specific stress-tensor Jmn under consideration. Accordingly, the met- 
ric 'perturbation' g = {R'^/z'^)h, which describes the shockwave, needs not be parametrically 

small. In fact, the normalization of this perturbation is related to that of the energy-momentum 
tensor on the boundary, via holographic renormalization [48,49]. Specifically, if h^'^\x~,x±) is 
the coefficient of the 2;^ term in the near-boundary (z 0) expansion of h(z,x~ ,x±), then 

r__(x-,xx) = £^^h(^\x-,x^) = ^/i(4)(x-,xx), (2.5) 

where the second equality follows after using the AdS/CFT correspondence to identify G5 = 
'kR^/2N^. As anticipated, T must scale like for the metric perturbation not to be para- 
metrically small. 

We shall now successively consider the two interesting cases. In this process, we shall recall 
some results from the literature, and we shall correct and extend the proposal in Refs. [12,30]. 

2.1 Type I metric: no source in the bulk 

In this case we need the non-trivial solutions to the homogeneous version of Eq. (2.4). Consider 



ffi'st the case where there is no dependence upon x±, corresponding to a nucleus which is uniform 
in impact parameter space (an infinite, uniform, wall). Then the solution reads 

97r2 

Hz,x-) = —z'T^^ix-), (2.6) 

c 

where the normalization has been fixed according to Eq. (p.5|). With this expression for h, the 
metric ( ^.3| ) is the original shockwave metric proposed by Janik and Peschanski [12]. Although 
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an exact solution to the (homogeneous) Einstein equations, this metric has nevertheless the 
drawback to have a singular point at z — > oo. We shall return to this issue later on. 



Consider now the homogeneous version of Eq. (2.4) corresponding to a generic profile in 
x±. It is then convenient to perform a Fourier transform to transverse momentum space, which 
yields the following equation 

[d^,-^d,-kl^h{z,x-,k^) = 0. (2.7) 
The general solution is expressed in terms of modified Bessel functions of second rank: 

h{z,x-,k^) = ^{zk^f[ci{x-,k^)K2{zk^) + C2{x-,k^)h{zk±)]. (2.8) 
The near-boundary expansion of this solution reads as follows (up to order z^) 

h{z,x-,k^) = ci{x-,k^)[l-^ + ^{-2ln{zk±/2) + 3/2-2jE) 

+ C2{x~,k±)^ +■■■ , (2.9) 
lb 

that is, the component proportional to l2{zk±) vanishes like z^ when z ^ 0, while that propor- 
tional to K2{zk±) approaches a non-zero value in this limit. Consider also the large z behavior: 
for zk_i_ ^ 1, K2{zk±) is exponentially decreasing, while l2izk±) is exponentially increasing. 
For what follows, it is interesting to keep in mind that the homogeneous (in the sense of no 
dependence upon x±) solution in Eq. ( |2.6|) corresponds to the limit k± ^ of the l2-piece of 
the general solution Eq. (|2.8|) . 

In the framework of AdS/CFT, the solutions which diverge as z — > oo are unacceptable and 



must be discarded. In view of that, it might look natural to enforce C2 = in Eq. (2^): this is 
the common strategy for computing correlation functions in the vacuum of = 4 SYM (here, 
the correlators of T^i,) [7-9], and this was also the proposal made in Ref. [30] for constructing 
the gravity dual of a nucleus. However, in what follows we shall argue that choosing C2 = 
in Eq. ( ^^ ) leads to a physically unacceptable picture for a nucleus. (In particular, this would 
also exclude the homogeneous shockwave ( |2.6D , which as alluded to above represents the limit 

— > of the piece of the solution proportional to C2.) A more sensible choice, which is 
physically motivated, is to take ci = and introduce a cutoff in the radial direction of AdSs at 
a distance z\ = 1/A, with A playing the role of an infrared cutoff in the boundary gauge theory. 
To motivate this proposal, let us first explain the difficulties with the original choice in Ref. [30] . 

As just mentioned, taking C2 = is the standard choice for computing the correlators of T^^u 
in AA = 4 SYM. In that case, the metric perturbation has a non-zero limit on the boundary, 

- = Sf,^5u-ci{x~ ,k±) , (2.10) 



hnu{x ,k±) = lim 

2-^0 



j^2 9fiiy{z, X , A; J 



which acts as an external source (not to be confused with the bulk source Jmn in Eq- ( |2.2| )), 
which couples to T^y in the boundary gauge theory. This external source induces a non-vanishing 
expectation value for T^^,, which is proportional to the source and can be read off Eq. ( ^ ) : 



647r2 



(T^.(x-, fc±)) = 5^-5,^ ^1 In tV ci(x-, k^) . (2.11) 
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Here, n is scale for ultraviolet renormalization in the gauge theory, as introduced by the removal 
of the logarithmic singularity at z — > manifest in Eq. ( |2.9| ). (The finite terms beyond the 
logarithm which are also visible in Eq. ( |2.9D have been absorbed in the definition of By 

taking a functional derivative in Eq. ( p. Ill ) with respect to h = ci, one obtains the only 

non-vanishing 2-point function of T^^, within the present kinematics (in momentum space and 
with trivial delta functions removed) : 

(r__(A^±)r++(-^±)) = ^kiln(^) . (2.12) 



647r2 ^ V^i, 

This is indeed the expected result^, with the specific /c^ -dependence in the r.h.s. reflecting the 
conformal symmetry of = 4 SYM (see, e.g., the discussion in [8]). 

Within this standard procedure, the external source h^i, plays no dynamical role, but rather 
acts as a device for generating vacuum correlations via functional differentiation. By contrast, 
in Ref. [30] it has been proposed to use a similar procedure with a suitable choice for the 
function ci{x~,k_i) in such a way to generate, via Eq. (2.11), an energy-momentum profile 



which physically would correspond to an ultrarelativistic nucleus. However, Ref. [30] showed no 
explicit proposal for such a function ci{x~ , k±), and in fact we shall now argue that there is no 
meaningful solution of this type. 

Our main objection to the proposal in Ref. [30] is of physical nature, but it also gets 
reflected in severe technical difficulties, which lead to paradoxes. Let us start with one of these 
paradoxes, before we explain the deeper origin of the problem. Returning to the simpler notation 
T = {T ), we notice that one must have T {k± = 0) > because 

r__(x^,/c_L = 0) = J d'^x±T__{x~,x±) (2.13) 

where T {x~,x±) is the positive semidefinite energy density in the transverse plane. For this 

to be consistent with Eq. ( ^.11[) one should have ci behaving like l//c^ when k± 0. But 
then the Fourier transform of h{z, x~ , kj_) back to the transverse coordinate space is ill defined, 
because of an infrared divergence at kj_ 0. On the other hand, there is clearly no similar 
problem if one chooses ci = and C2 7^ in Eq. (|2.8D , since then one can take C2 ~ 1/fc^ as 
A;_L ^ without generating infrared problems, as obvious from the expansion (|2.9| ). 

That this is the only meaningful choice, it is also suggested by the behavior (|2.9[) of the metric 
near z = together with its physical interpretation according to the AdS/CFT dictionary. Via 
the UV/IR correspondence [26,31-35], the inverse l/z of the radial distance in AdSs is mapped 
onto the transverse momenta (or virtualities) of the quantum fluctuations in the boundary 
gauge theory. Hence an expression like Eq. ( |2.8| ) for the metric perturbation should be viewed 
as encoding information about two types of transverse momenta: the momentum k± which via 
the Fourier mode C2{x~,k±) determines the profile of the nucleus in impact parameter space, 
and the momentum ~ 1/z which refers to the quantum modes that we would like to include 
in our description of the nucleus and of its high-energy interactions. 



^In the corresponding expressions in the Uterature, k'^ is generaUy replaced by the invariant 4-momentum 
squared = k^— 2k^k~ . Note however that for the present set-up we have k~ = since there is no dependence 
upon . 
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On physical grounds, in particular in view of our experience with QCD, we expect the bulk 
of the modes in a hadron wavefunction to be concentrated at 'soft' momenta, of the order of 
some infrared cutoff A, whereas at much larger momenta p± ^ A the distribution should rapidly 
decrease, typically according to a power law The exponent A is equal to 2 for the 

partonic tail produced via bremsstrahlung in QCD at weak coupling, but it is equal to 4 or 
larger in all the situations where a gravity dual has been identified for a hadronic system at 
strong coupling. For instance, A = 4 in the case of the AdSs black-brane geometry dual to 
the Af = 4 SYM plasma [11], and also for the uniform shockwave metric in Eq. (|2.6D [12], or 
for the metric induced by a pointlike source in the bulk of AdS^ [17] (see Eq. (|23^ ) below). 
Furthermore, for a normalizable dilaton state in AdSs [8,9], which is dual to a 'glueball', one 
has A = 2 + \/ m?R^ + 4 > 4, where the dilaton mass m is proportional to the scale A which 
enters via the radial cutoff at z ~ 1/A [23] . Such a large value for A at strong coupling can be 
understood as the absence of partons with high momenta [33-36] , an interpretation to which we 



shall return in Sects. 3 and 4. Clearly, A = 4 is also the prediction of Eq. (2.8) provided one 



takes ci = (cf. Eq. ( |2.9D ). On the other hand, with C2 = 0, Eq. ( |2.8| ) predicts a dominant 
behavior near z = which is independent of z, that is, a flat momentum spectrum (A = 0) in 
the high momentum regime, which physically makes no sense. 

The previous examples also show that the behavior in of the metric perturbation holds 
only for sufficiently small values of z, corresponding to the ultraviolet behavior of the dual gauge 
theory. But this growth is cut off at some larger value of z, whose inverse plays the role of an 
infrared cutoff in the gauge theory. Such a cutoff is necessary to introduce the analog of matter 
(plasma or hadronic bound state) in the otherwise conformal SYM field theory. For instance, 
this scale is provided by the black hole horizon at z ~ 1/T for the AdSs-Schwarzschild metric, 



by the radial position of the source in the bulk for the shockwave metric in Ref. [17] (see Sect. 2.2 



below), and by an explicit cutoff at large values of z in the construction of the glueball bound 
state [22,23,25]. 



At a first sight, the uniform shockwave metric ( |2.6| ) seems not to involve any such a scale, 
but this is only illusory: written as it stands, the metric perturbation in Eq. ( |2.6D blows up 
at z — > oo and this divergence is a genuine space-time singularity, and not just a coordinate 
singularity. As mentioned in the Introduction, this singularity is quite elusive, as it does not 
show up in any of the curvature scalars computed from this metric. Rather, it can be identified 
via a more general criterion for space-time singularities, known as (7-incompleteness (with g 
standing for "geodesic") [29]. Namely, in an acceptable space-time, which is (7-complete, any 
time-like or null-like geodesic can be extended up to arbitrary values in their affine parameters 
(like the proper time). The contrary would imply the existence of observers whose history ends 
or begin at a finite proper time. On the other hand, for the Janik-Peschanski metric one can 
construct time-like geodesies which, starting at a generic point zq, reach infinity after a finite 
value of the proper time. (See [27] for an analogous construction in the Kaigodorov space- 
time, which is the 4-dimensional version of the Janik-Peschanski metric.) Hence, this metric is 
^-incomplete, and thus unacceptable. One can effectively 'hide' the space-time singularity at 
z ^ 00 by introducing a radial cutoff at z = 1/A, which for that purpose plays the same role 
as the black hole horizon within the AdSs black-brane metric. But once this is done, it is no 
more disturbing to use a general metric perturbation cx z'^l2{zk±), which has the right behavior 



- 11 - 



at small z and allows for a generic transverse inhomogeneity, instead of the uniform shockwave 
in Eq. pIBI). 



To summarize, our proposal for a gravity dual to a large nucleus is given by the asymptoti- 
cally AdSs metric in Eq. ( p. 3D which applies for z < za = 1/A and where the function h is given, 
in transverse momentum space, by Eq. (|2.8| ) with ci = and C2 related to the energy-momentum 
tensor of the nucleus that we would like to describe via 

r__(x-,A;^) = ^fcic2(x-,A;x). (2.14) 

The example of the 'plasma slice' previously discussed suggests that a physically reasonable 
choice for T would be 

T__(x-,xx) = 7'A^cA'/(x-,xxA), (2.15) 

where the dimensionless function / describes the shape of the 'nucleus' in longitudinal and 
transverse directions. The transverse inhomogeneity is controlled by the 'soft' scale A, so like in 
QCD. The longitudinal support is concentrated at < L/7 with L the width of the 'nucleus' 
in its rest frame. In the high energy limit, we can replace this by a (5-function in x~ : 

f{x-,x^A) = f{x^A)-6{x-), A2 [ dV/(x±A) = 1, (2.16) 

7 J 

where the normalization of the new function f{x±A) has been chosen for convenience. Some 
reasonable choices for /, inspired by our experience with QCD, would be a Gaussian, or a 
Woods-Saxon profile which falls exponentially for distances far away from the center, that is for 
xj_A ^ 1. These choices lead to the following model for the metric perturbation: 

h{z,x-,x^) = m^^^LA^zH{x-) I e^''^^^f{k^/A)h{zk^), (2.17) 

where f{k_i_/A) is the Fourier transform of A^/(xj_A) (so it is dimensionless), and it has support 
^-t k± A. A radial cutoff at z ~ za is implicit. Then, clearly, the argument of I2 can never 
become large. 

Our prescription for cutting off AdS at z = z\ is merely a model, and any calculation which 
is sensitive to large values of z (so like the scattering between two Shockwaves) will be strongly 
sensitive to the details of this model — the value of A and the specific procedure used for its 
implementation. Fortunately, there are also interesting phenomena, so like the deep inelastic 
scattering to be considered in Sects. 3 and 4, which are controlled by the 'hard' {p± ^ A) part 
of the spectrum — in the AdS framework, by the behavior of the metric near the boundary at 
z = — and thus are completely insensitive to the model used to cutoff AdSs at large z. 

2.2 Type II metric: pointlike source in the bulk 

We now briefly describe the shockwave proposal in Ref. [17], which involves a source in the bulk. 
We shall take this source to be composed of point-like 'particles' moving together* along 



*Ref. [17] mentioned only a single such a particle, but here we shall consider a collection of A'^^ of them, in 
order to achieve a metric perturbation of order one. The alternative possibility, which would be to take a single 
particle but with ultrahigh energy p"*" ~ , would be inconsistent with the use of the supergravity approximation 
for any collision involving that 'particle' [33,36]. 
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the trajectory defined by: x = 0, x_l = 0, and z = z^, with a positive constant. Then the 
only non-zero component of Jmn is J , and is given by [17] 

J__ = p+N^ ^ 6{x-)6'^^\x^)6{z - z,) , (2.18) 

with p'^ the light -cone longitudinal momentum of a single 'particle'. The solution to Eq. ( p.4| ) 
corresponding to this current reads 

h(z,x-,x±) = -^^^2i^i(3,5/2,5,-l/(7)<5(x-) (2.19) 
16 q-^ 

where 

_ x\ + {z- z^f 



4:ZZ^ 



(2.20) 



is the so called 'chordal distance'^ between the point {z,x±) where we measure the metric and 
the position {z^,x± = 0) of the pointlike source in the bulk, and the hyper geometric function 
takes a rather simple form: 



3 1 + 8(7(1 +q)- Vg(l + <7)(1 + 2q) 



2^1(3,5/2,5,-1/(7) = I6q' ^. (2.21) 

As anticipated, the energy E of the particle must scale like in order for the metric pertur- 



bation Eq. ( 2.191) to be non-negligible. Note that g 3> 1 for both small (z ^ z^,) and large 



[z ^ z^) values of z. Using 2F1 ~ 1 — 3/2q 1 when g » 1, we deduce 



h{z,x-,x^) ^ ^^^^^j^6{x-) for z«z., (2.22) 

[X I -r Z^ ) 



and respectively 



47rp+(zz*)'^ 

{xl + Z2)S 



h{z,x ,x±) ~ -j--^ — ■ — ^S{x ) for 2; » z,, , (2.23) 



so, in particular, the metric perturbation dies away as 1/z^ when z — > 00. 

The associated energy-momentum tensor in the boundary gauge theory follows from Eq. ( |2.5| 



0+ /v2 2z'^ 
vr (x^ + Z;^)^ 

This is essentially uniform so long as xj_ <^ z* but it decreases like 1/x^ for xj_ ^ z*. We 
see that 1/z* plays the same role as the 'soft' momentum scale A introduced previously, in the 
sense of fixing the scale for transverse inhomogeneity in the nucleus. But unlike in that previous 
case, now there is no need to explicitly cut off the radial dimension of AdSs, since the metric 



perturbation Eq. ( 2.19 ) dies away, like 1/z^, when z — > 00. 



^The chordal distance is the 50(3, l)-invariant distance, i.e. the analog of the radial distance, for the hyper- 
bolic space H3 spanned by the coordinates {z,x±). 
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3. DIS off the Shockwave: the dilaton 



In this and the next coming section, we shall consider the supergravity problem dual to the deep 
inelastic scattering (DIS) between a 'hard' (i.e., highly virtual) external current and a nucleus 
in the M = A SYM theory. The 'nucleus' should be thought off as a slice of a plasma and it will 
be described as a shockwave, as already explained. Its construction requires an infrared cutoff 
A, but the details of this cutoff will be unimportant for the hard process at hand. The 'current' 
generally associated with DIS is a virtual photon with space-like virtuality which couples to the 
electromagnetic current of the nucleus. Within = 4 SYM, this 'electromagnetic current' can 
be represented by the conserved T^-current carried by (adjoint) fermionic and scalar fields in 
the Lagrangian. The DIS of such an 7?.-current will be addressed in Sect. But before doing 
that, it is preferable to introduce the formalism in the simpler context of the scalar 'current' 
J = {l/A)F^^Fa^ (a is the SU(A''c) color index). This operator is interesting in itself, since it 
couples to the gluons inside the nuclear wavefunction, and hence is a direct probe of the gluon 
distribution. The 'bulk' AdS field dual to this operator is the massless dilaton field, which obeys 
the simplest equation of motion in supergravity: the Klein-Gordon equation in the relevant 
(here, the shockwave) metric. For simplicity, we shall refer to the projectile as the 'dilaton' 
although, strictly speaking, it is the operator J, and not the dilaton, which undergoes DIS in 
the boundary gauge theory. 

3.1 Formalism and structure function 

As before, we shall work in the infinite momentum frame of the 'target' (the nucleus), which 
therefore will be taken to be Lorentz contracted to a (5-function at x~ = 0. The 'projectile' (the 
dilaton) propagates in the negative direction (a 'left mover'), with space-like momentum q'^. 
In light-cone coordinates, we have > and < 0, and we take q_L = for convenience; 
hence, the virtuality of the projectile reads = q'^q^ = —2q'^q~' > 0. As usual with DIS, it 
is convenient to express the total cross-section (or the 'structure function') in terms of and 
the Bjorken-j; variable, defined as 

X = = = (3 1) 

^ -2q-p 2q~p+ 2q--fA ^ ' ' 

where = d^'^p'^ is the momentum of a typical constituent of the target: p"*" = 7 A with A 
the characteristic momentum scale in the target rest frame (e.g., A = T for a slice of a finite- 
temperature plasma) and 7 1 the nuclear boost factor. The kinematic conditions of interest 
are such that Q'^ > A^ and x <C 1. 

Via the optical theorem, the DIS structure function F{x,Q'^) is expressed as the imaginary 
part of the forward scattering amplitude for Compton scattering: 

F{x,Q^) =lmU{x,Q^), (3.2) 

with the 'polarization function' 11(2;, Q^) usually written as 

n(x,Q2) I dVe-^''-(P|T{J(x)J(y)}|P), (3.3) 

where r'^ = x^—y^. However, as above indicated, this writing is appropriate only when the target 
is in a state | P) with given 4-momentum P^, which is translationally-invariant. Here, however. 
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the target is represented as a shockwave which is locahzed in space and inhomogeneous in both 
longitudinal (x~) and transverse {x±) directions. Accordingly, the current-current correlator 
computed in this shockwave background depends not only upon the relative separation r'^, 
but also upon the central coordinates b~ and b± (defined as = (x^ + y^)/2). To obtain a 
structure function which depends only upon the kinematical variables x and Q^, we shall follow 
the prescription in Ref. [50] and average out the central coordinates. The longitudinal extent of 
the target is of order L/7 (recall that L denotes its width in its rest frame) whereas the scale 
for inhomogeneity in the transverse plane is set by A. We shall thus replace Eq. (^]^) by (the 
precise normalization is irrelevant at this point) 

n(x, Q2) = (7/L)A2 y dV db- d'b^ e"i ' iG(xo) {[J{b + r/2), Jib - r/2)]) , (3.4) 

where the expectation value is now computed according to the AdS/CFT correspondence (see 
below) and we have replaced for convenience the time-ordered correlator by the retarded one. 
Both types of correlators yield the same imaginary part, but the retarded one is more easy to 
evaluate in the context of AdS/CFT; see e.g. [3]. Note that J cx {F^y)"^ has mass dimension 4, 
so n will have dimension 4. 



Let n(x,y) denote the (retarded) 2-point function of J which enters Eq. (3.4). According 
to the AdS/CFT correspondence, this is formally obtained as 

n(x,w) = ^ , ,^^f, , X (3.5) 
where Sd is the action of the dilaton field, that is, 

5 = -TIT^ / A^xdz^gg^^dn<\>d^<^, (3.6) 



167r2fl3 

evaluated with the solution ^{z^ x) to the classical equations of motion, 

5™(^/^9""9„'A) = 0, (3.7) 
obeying the following boundary condition at 2; = : 

lim</)(z,x) = 4>h{x) . (3-8) 

On physical grounds, we need the boundary field 0b (x) to be the plane-wave 

0b(x) = e-^(''"-++''+-") 0b (3.9) 

(0b is simply a number), but at intermediate steps we shall consider a generic function (pbix), 
for the purposes of functional differentiation. Using the equations of motion, one can perform 
the integral over z in Eq. (|3.6| ) and obtain 

j '^'^^ -^Hz,x)d^<p{z,x)\^^^ (3.10) 



ri 



167r2 

Note that the boundary condition (^.8]) is not enough to fully specify the classical solution. 



since Eq. (3.7) is a second order differential equation in z. Furthermore, the above procedure 
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cannot reproduce the imaginary part of the retarded correlator, as obvious from the fact that 
the 2-point function generated via Eq. ( p.5D is symmetric in its arguments. We shah return to 
these problems later on. 



Eq. (3.7) is written in the shockwave metric ( |2.3| ), which is independent of the light -cone 
time x"*". Thus, the equations are homogeneous in x"*", meaning that the variable q~ introduced 
by the projectile is conserved by the dynamics. It is therefore sufficient to restrict ourselves to 
the respective Fourier mode, by writing 

(f){z,x~^ ,x~ ,x±) = e~'^ ^'^ (l){z,x~ ,x±). (3.11) 

(The function (j){z,x~ ,x±) in the r.h.s. depends upon , but this dependence is kept implicit.) 
Then the equation of motion is more explicitly written as (with the notation x = {x~ ,x±)) 

- + 2iq-d^ + Vi) ^{z, x) = -{q-fh{z, x) 0(z, x). (3.12) 

We have separated in the right hand side the interaction piece, which describes the scattering 
between the dilaton and the shockwave. We shall correspondingly decompose the total field cp 
into its 'incoming', or 'vacuum', piece (po and the 'scattering' piece cps'- 

(piz,x) = (j)o{z,x) + (j)siz,x) , (3.13) 

We construct (j)o in such a way to satisfy the boundary condition ( |3.8D , hence 

lim(j)Q{x, z) = (j)b{x), lira (j)s{x, z) = 0. (3.14) 



The vacuum problem is formally similar to the source-free version of Eq. (2.4) that we have 
already solved: the equation is homogeneous in space and time, so the solution corresponding 
to the boundary field ( pl9| ) is of the form (l)o{z,x) = e~*^^^' (l)o{z). The function (poiz) obeys 
the equation obtained by replacing k'j_ — > —2q'^q~ = > in Eq. ( ^ ) with the boundary 
condition 0o(z ^0) = 4>b- Clearly the unique acceptable solution is (compare to Eq. (2.8)) 

M^) = l{Qzf'i^2{Qz)^b- (3.15) 

This solution, together with Eqs. ( |3.5| ) and ( p. 10 ), yields the vacuum component no(Q^) of the 
polarization function in a form entirely similar to Eq. ( 2.12| ) : 

no(Q') ^^ J d^xe-«-(0|T{J(x)J(0)}|0) = ^ In . (3.16) 

For what follows, it is useful to rewrite Eq. ( p.l5| ) in a more general way, which features the 
vacuum boundary-to-bulk propagator Dq{z,x — y) : 

4>q{z,x) = j d^yDo{z,x-y)My), Do{z,k) = ^{KzfK2{Kz), (3.17) 

where the momentum-space version of Dq was written for the space-like kinematics of interest 
here: = —2k~^q~ + k'j_ > 0. Via manipulations to be described in the next section, we will 
construct the scattering piece (ps in a similar form: 

0,(z,x) = / d'^yD{z,x,y)My), (3-18) 
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where D{z,x,y) is the boundary-to-bulk propagator in the shockwave metric ( |2.3| ) and is in- 
homogeneous in the spatial coordinates. Inserting this into Eq. ( p. 10 ) we obtain the scattering 
piece of the classical action as 



Sd — So 
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d^x / (l^y(j)b{x)(j)b{y)-^dzD{x,y,z 



(3.19) 



When acting on ( 3.19 ) with the functional derivatives in (p.5|), we obtain two terms: 



6(t)b{x) 6(l)biy) 



d^x / dc^y (l)b{x) My)^ 9zD{x,y, z) 



2 = 



{D{x, y, z) + D{y, x, z)) 



(3.20) 



which together would provide a symmetric and real expression for the current-current correlator 
n(a;, y). Clearly, this is not the physical result that we are interested in. To recover the imaginary 
part of the retarded 2-point function, and hence the structure function -F(x, Q^), we follow the 
prescription in Refs. [51,52] and drop the second term in Eq. ( |3.20| ), while multiplying the 
contribution of the first term by two: 

^R{x,y) = — d;,D{x,y,z) 



167r2 



(3.21) 



Note that, as it will become explicit later on, the propagator D{x,y,z) is retarded with respect 
to the variable x~ , which plays the role of the light-cone time for the dilaton field (a left-mover). 



Moreover, the solution built with this propagator, cf. Eq. ( 3.18 ), is such that for large times x~ 
it represents a purely 'infalling' wave, i.e. a wave which with increasing x~ propagates towards 
larger values of the radial dimension z. Thus our above prescription for keeping D{x, y, z) while 
discarding D{y,x,z) is indeed the same as the prescription for keeping the infalling solution 
alone, as originally formulated in Refs. [51,52]. 



Given Ilji{x,y), the polarization function (3.4) is finally computed as 



n(x,Q2) = (7/L)A^ / dx-dy- d'^x^d^y^e 



Jq+{x -y ) 



^R{x,y) 



(3.22) 



where as compared to Eq. (3.4) we have removed the Fourier transform over x~^ (since this is 
automatically performed by working in the (/"-representation) and rewritten the measure as, 
e.g., dr^dft^ = dx^dy^. One can check on the previous equations that n^(a;,y) has mass 
dimension 7, hence Il.{x.,Q'^) has dimension 4, as it should. 

In Sect. 3.3, we shall find that the structure function F{x,Q'^) = lmIl{x,Q'^) obeys the 
following sum-rule 



hm ^ [\xxF{x,Q^) = ^N^A' 



(3.23) 



which is recognized as the statement of energy-momentum conservation: this particular integral 
of F{x,Q'^) singles out the coefficient in front of T^i, (the energy-momentum tensor operator) 
in the operator product expansion (OPE) of the current-current correlator. The quantity in the 
r.h.s. of Eq. ( |3.23 ) can be recognized as the nuclear energy density in its rest frame (up to a 
normalization factor). 

^"it is now understood that the functional derivatives are defined with respect to the three-dimensional field 
(f>b{x' ,x±) and for a fixed value of q~ . 
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3.2 Eikonal scattering 

In what follows we shall construct the scattering field (ps by iterating the interaction piece in 
the r.h.s. of Eq. ( |3.12| ). This amounts to resumming graviton exchanges between the target and 



the projectile to all orders. To that aim, we need the integral version of Eq. (3.12), that is, 

/di 
-p dy'<i^y± G{z, z; x~ -y~,x±- y±) 

X [-iq-f]hiz,y-,y^)(l)iz,y',y^), (3.24) 
where G{z, z; x~ — y~ ,x± — y±) is the bulk-to-bulk propagator obeying 



- - 5. + 2iq-d. + Vi j G{z, z; x--y-,x^- y^) = zH{z - z)5{x- - y-)6^^^ {x^ - y^) , 

(3.25) 

to be explicitly constructed in Appendix |B|. (The euclidean version of this propagator is well- 
known in the literature, and will be recovered in Appendix but here we rather need its real- 
time version in light -cone coordinates and mixed Fourier representation.) Given the boundary 



conditions Eq. (3.14), it is clear that G must vanish at the boundary. Several expressions for 



this propagator will be useful in what follows. For instance, the following one 

^2i.. /-oo 



Giz,z;x -y ,x^-y^) = / (^ « 



dujuj 







X z2j2((jz)i2j2(cji)e ^1" , (3.26) 

makes it clear that the propagator vanishes as z — > and/or i ^ 0, and also that it is retarded 
with respect to the projectile time variable x" — y^ ■ Furthermore, by using the completeness 
relation in Eq. (|B.3|) for the Bessel functions, one can show that 



Giz, z; X- -y- ^ +0, x^-y^) = --^ z^6{z - z)5^^^ {x± - yx) • (3.27) 

2q 



We now proceed to formally solve Eq. ( |3.24| ) via iterations. The first iteration gives 

di 

X h{z,x~ ,x±)(l)o{z,x~ ,x±), (3.28) 



')i{z,x ,x±) = —{q )^ / ~ d^x^ G{z, z;x —x ,x± — x±) 



which we shall compactly rewrite as 

Mx) = -{q'f J dXG{x,x)h{x)Mx). (3.29) 

The second iteration gives 

and so on. An important simplification occurs in the limit where the shockwave h is treated as a 
(5-function in x~; then, the perturbative series produced by iterations exponentiates and yields 
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an eikonal phase — a natural result at high energy. Specifically, let us write h{z,x ,x±) 
6{x-)h{z,x_L). Then from the factor h{X)G{X, X)h{X) in Eq. ( |3^ ) we obtain 



6{x )Q{x —X )6{x )= -S{x —X )6{x ) 



which after also using Eq. ( |3.27| ) gives a factor 



' -S{X-X)6{x 



2q- 2 



Thus Eq. ( |3.30| ) reduces to 



2q- 



AX i X) 5{x-) h\X) MX). 



(3.31) 



(3.32) 



(3.33) 



This procedure is easily generalized to higher orders: for the k-th. iteration we deduce 

k-l 



1 



{q-f]^h\X) 



2q' 



(3.34) 



where the factorial is again generated via the product of theta functions in the propagators 
and delta functions in the metric field h. Then the sum exponentiates, as anticipated, with the 
following final result: 



[z,x ,X_l) 



/dz 
— d^x± G{z, z; x~ ,x± — x±) T{z, x±) 4>o{z, 0, x±) 



(3.35) 



where we have defined the scattering amplitude corresponding to a radial penetration z for the 
dilaton and in the eikonal approximation as 



iT(z, b I 



1 — exp 



h{z,b±) 



(3.36) 



An interpretation of this quantity in the boundary gauge theory will be shortly given, following 
Eq. ( |OT|) . 

Using the expression of (ps in Eq. ( ^) it is now straightforward to compute the DIS 
structure function according to Eqs. ( |3.2| ), ( 3.21 ) and ( 3.221 ). First we see that the bulk-to- 
boundary propagator defined in equation ( |3.18D is given by (cf. Eq. ( 3.17]) ) 

D{x,y,z) = -2q' j d^b^ j ^ G{z, z; x~ ,xj_ - b±)T{z,bj_)Do{z, -y' ,bj_ - y±). (3.37) 

Clearly the integrals over x~ and y^ in Eq. ( p. 22 ) select the Fourier components with = q^ 
in both G and Dq. So we are left with 



n(x,Q2)=-^^2g~y dVdV I J ^^d,G{z,z;q+,x^-b±) 



z=0 



xT{z,b^)Doiz,q+,b^-y^). (3.38) 
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Using the transverse momentum space form ( [B.5[ ) for the Green's function, we see that the z 
piece in the above reads 

2 



-{KzY^^iKz) 



(3.39) 



2=0 



where = ^1 ~ = A;^ + Q^. By also using Eq. ( |3.17| ) for Dqi we deduce 



n(x,Q2 



2g- 



327r2L 
{Kzf¥.2(Kz)T{z,h 



d^rc^d^y^ / d^6^ / ^ 



dz /• d^fcx 
(2^ 



Jfcx-(a:x-6±) 



(Qz)2K2(Qz), 



(3.40) 



(27r)2 

with = + This expression simplifies drastically after performing the transverse 
coordinate integrals, which set A:j_ = = 0. Recalling the definition Eq. (^]^) we finally arrive 
at 

n(x,Q2) = y dzzK2(Qz) y dVr(z,6x), (3.41) 

from which the structure function -F(x, Q^) finally obtained by taking the imaginary part. 

This result has a natural interpretation: K^Q^) plays the role of the wavefunction squared 
for a partonic fiuctuation of the 'current' J = (l/4)i^2 which according to the UV/IR corre- 
spondence has a transverse size r_|_ ~ z. The modified Bessel function effectively restricts r_|_ to 
values ~ z < 1/Q, as expected from the uncertainty principle. Furthermore, 



a{z,x) = 2 j d%_LlmT{z,b±) 



(3.42) 



is the total cross-section for the scattering between this partonic fluctuation and the 'nucleus'. 

Let us finally verify that the use of the (5-function approximation for the metric perturbation, 
and hence the eikonal approximation, are indeed justified for the problem at hand. By inspection 
of the previous manipulations, it is clear that the only place where the assumption that h oc S{x~) 
has played a role was in using 



2q- 



< 1 



(3.43) 



in order to simplify the last exponential in Eq. ( 3.26| ) and thus replace the intermediate propa- 
gators in the iterative series by (5-functions, cf. Eq. ( |3.27| ). In reality x~ is, of course, not strictly 
zero, but rather of order L/7, which is the longitudinal width of the target in the infinite mo- 
mentum frame. Since k± can never be too large (this is set by the transverse inhomogeneity in 
the target, so ~ A), while uj is typically of order 1/z (as fixed by the oscillatory behavior of 
the Bessel functions in Eq. ( p. 26 )), the above inequality amounts to 

z' » -— = . 3.44 

2q 7 Q2 

Note that, in this argument, z is the radial distance at some generic scattering point, e.g. z 



or z in Eq. ( 3.30 ), and not the argument z of in Eq. ( 3.35 ), which approaches zero when 
computing the structure function, cf. Eq. ( |3.39| ). As we shall shortly check, the integral over 

^ l/Q^ ; hence, the condition in Eq. (3.44) is well satisfied 



in Eq. ( 3.41 ) is dominated by 



whenever xAL <C 1, which is indeed the interesting situation (since, typically, L ~ 1/A). 
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3.3 Prom single scattering to the saturation momentum 

To further compute the DIS cross-section according to Eq. ( |3.41| ) we need to specify the metric 
perturbation h{z,b±) which enters the scattering amplitude ( |3.36D . The crucial point for what 
follows is that, as manifest on Eq. ( |3.41 ), the DIS process is only sensitive to relatively small 



values of z, such that 2; < 1/Q <C 1/A. In view of this and of the discussion in Sect. it is clear 
that, for this purpose, it is enough to retain the dominant behavior of h{z, b±^) near the boundary 

{z ^ 0), which for any acceptable shockwave metric scales like z'^T , cf. Eq. (p.5|). Thus our 

formalism makes it clear that the DIS process at strong coupling involves the scattering off the 

nuclear energy density T On the supergravity side, this scattering involves multiple graviton 

exchanges (here resummed in the eikonal approximation), with each factor of corresponding 
to one exchanged graviton. Hence, from now on we shall simply take 

J" 

h(x-,b^,z) = h(b±, z)S(x-) = 27r2— - , (3.45) 

which is the approximate solution to the Einstein equation for the metric component h so long 
as zA <C 1. One can cover both type of Shockwaves described in Sect. |2| by taking 

T^-{x-,b±) = k^Ef{b^K)5{x') = f^-{bi_)5{x") , (3.46) 

where E is the total energy of the nucleus: for the type~I shockwave in Eqs. ( p.l5D -( ^.16| ), this 



reads E = ■jLA N^, whereas for the type~II one in Eq. ( 2.24 ) one has E = p and we identify 



= 1/A. 

Let us now consider the single scattering approximation which amounts to expanding 
T{z, b±) in Eq. ( 3.36| ) to lowest order. By also making use of Eq. ( 3.45| ) the amplitude reads 



rW(.,6x) = g-v^2^^.^ (3.47) 

which is purely real, reflecting the fact that the single graviton exchange generates no imaginary 
part (one cannot cut through a single graviton propagator). Note that the single graviton 
exchange is proportional to the energy density in the projectile and the target, and is of order 
zero in A'^^, because the target contains degrees of freedom per unit volume (otherwise it 
would be suppressed like l/N^). Via the UV/IR correspondence z <^ r±, Eq. ( |3.47| ) implies that 
a small dipole fluctuation of the projectile, so like a gluon-gluon pair, with transverse size r± 
interacts with the target with an amplitude T{r±) ~ r^, which vanishes much faster at small 
r± than at weak coupling^^ (where one rather has T ~ [6]). As we shall later argue, this 
difference is due to the fact that, at strong coupling, there are no high-Q^ partons in the target 
wavefunction, that the dipole could scatter off. 

In this single-scattering approximation the polarization function simplifies to 

n«(-,Q^) = ^^/ d.dV.Ki(Qz) ''^^/^^^ z\ (3.48) 



^^This would not be true for the shockwave metric selected in Ref. [30], and which is obtained by taking C2 = 
and ci / in Eq. In that case, the DIS cross-section would involve only the first, constant, term in the 



small-z expansion (2.9), which via the UV/IR correspondence would imply that the scattering amplitude remains 
constant as rj_ ^ 0. 
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which is purely real and hence describes only elastic scattering. The integrals in Eq. ( 3.4g| ) can 
be exactly perform and they yield the same final result for both types of metric (for the type-I 
metric, the integral over b± is performed according to Eq. ( 2.16P ), which reads 



10x2 



n«(x,Q2)= (3.49) 



This is suppressed by a power of A^/Q^ with respect to the corresponding vacuum result in 
Eq. ( p^ and thus it is recognized as a leading-twist effect. The fact that this is purely real, 
i.e. the structure function F{x,Q'^) vanishes in the leading-twist approximation, means that 
there are no point -like constituents in the nuclear wavefunction at strong coupling, in agreement 
with previous results in [33,34,36]. The l/x^ rise in Eq. ( 3.49| ), which on the supergravity side 
is clearly associated with the graviton exchange, can be also understood in the original gauge 
theory, as we explain now: in the OPE of the current -current correlator, valid at high ., 
the leading-twist operators with spin n should produce contributions which behave like l/x". 
At strong coupling, one expects all such operators to acquire large and negative anomalous 
dimensions [53-55], of 0{\^^^), with the exception of the spin 2 energy-momentum tensor 
operator, which is protected by symmetry. Hence, the only leading-twist contribution which is 
expected to survive at strong coupling is the one which behaves like 1/x^, in agreement^^ -^^ith 
Eq. (|3l9|) . 

Although it does not contribute directly to the structure function, the single scattering 



contribution in Eq. ( 3.49 ) allows us to derive the sum-rule ( p. 23)) . Specifically, let us introduce 
the variable u = l/x and extend the polarization function n(z/, Q^) to complex values for the 
variable v. Then, this function is expected to be analytic everywhere in the complex v plane 
except for two branch cuts along the real axis, at > 1 and u < —\ (the physical region for 
DIS and, respectively, the process related to DIS by crossing symmetry; see, e.g., [47]). Then, 
via contour integration in the complex v plane, one can relate the behavior of this function near 
= 0, where the OPE applies, to integrals ('moments') of the structure function F = Imll 
along the branch cuts. In the present context at strong-coupling, there is only one leading-twist 
contribution to Il{v,Q'^) — the single-scattering piece Ii^^\v,Q'^) oc cf. Eq. (|1|) — , 

meaning that there is only one moment of the structure function which survives in the high 
Q2 limit^'^. Clearly, this moment is proportional to the contribution of the energy-momentum 
tensor (the only protected leading-twist operator) to IT. This contribution is isolated as follows: 

— — = i ^ ~ 2 / ^ = - / dxxF(x,Q2 , 3.50) 

where the contour in the first integral is a small circle surrounding the origin which is then 
distorted in the complex plane in such a way to wrap around the two branch cuts which give 
equal contributions. The approximate sign in the second equality refers to the leading-twist 
approximation. As anticipated, this procedure reproduces Eq. ( 3.23| ). 

^^Via simple dimensional arguments, one can convince oneself that the contribution of Tfj,y ~ T--5^-Su- to 
the OPE of 11(2:, Q^) can be written as the product of /x^ times a quantity which represents the average energy 
density per unit transverse area in the nuclear target. This quantity is of order N^A'^, so the overall contribution 
has indeed the parametric form shown in Eq. (3.49). 

^''it is understood that one has to multiply the moment by the ratio A^/Q^ before taking the limit 00, 
so Uke in Eq. (^.23|). 
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The dominant contribution to the DIS structure function at high is of higher-twist order 
and comes from the two graviton exchange. This is easily estimated as 



F^'\xM^) = ^^^^ / d.dV.K^(Q.) 



2 



167r N'^K^L , 
K — 3.51 

7 x-^ 



with K = ttA^ J d'^t>±f'^; in particular, k = 1/8 for a type-I shockwave with the exponential profile 
/ = (l/27r) exp(— 6^A) and k = 4/5 for the type-II shockwave. Since generated by cutting in 
between two graviton exchanges, this contribution to the structure function can be characterized 
as being diffractive. Similarly, by expanding the eikonal phase in Eq. ( |3.36| ) to higher orders, 
one generates diffractive contributions to F{x,Q'^) of successively higher twist order and which 
increase faster and faster with decreasing x (as a power 1/x"^^ for the term corresponding to 
n graviton exchange). If one tries to use any single such a diffractive contribution within the 
sum-rule ( p.23| ) , the ensuing integral is ill-defined at its lower limit x ^ 0. This refiects the fact 
that the expansion of the eikonal amplitude makes sense only so long as the exponent is much 
smaller than one. Thus, clearly, one expects a change of regime when this exponent, which is 
the same as the single scattering amplitude, becomes of 0{1). 

To study this change of regime, let us use Eq. (|3.46D to rewrite 

(1) _ 2 EA^f{b±A) 4 A3l/(6^A) 



where we have also used the fact that z ~ 1/Q together with Eq. ( |3.1| ). (This estimate applies 
to both types of Shockwaves with the understanding that L ~ 1/A ~ for the type-II metric.) 
The condition that T^^) be of 0(1) can be solved for at fixed x and b±, thus defining the 
saturation momentum : 

Ql{xM) = ^^f{b±A), (3.53) 

where the factor 7r^/2 has been introduced for later convenience. Note that the transverse 
profile f(b±A) of the nuclear energy density directly transmits to the saturation momentum. 
In particular, for a shockwave which is homogeneous in the transverse impact parameter space 
one has = -KLh? /2x, in agreement with Refs. [36,40], but in disagreement with Ref. [39], 
where a rather surprising result for Qg was reported, which becomes independent of x at high 
energy. As it should be clear from the above analysis, and also from the corresponding ones in 
Refs. [36,40], the 1/x rise of at small x refiects the respective behavior of the one-graviton- 
exchange scattering amplitude, and thus it seems unavoidable within this gauge/gravity duality 
context, where the high energy scattering always amounts to graviton exchanges. 

Note finally that T'^^\z ~ 1/Q, 6_l) ~ Ql{x, b±)/Q'^, which reaches its maximal value at the 
center of the 'nucleus' (6j_ = 0). This makes it clear that the 'twist' expansion of the eikonal 
phase is an expansion in powers of Q^(x,0)/(5^ and is appropriate in the high-Q^ regime at 

^ Q'^{x,0). The opposite regime, at < Q'^{x,0), or equivalently x < Xs{Q) where, 

7^2 A3 r 

XsiQ) = ^g^/(0)> (3-54) 
will be studied in the next subsection. 
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3.4 Structure function at saturation 

Having identified the saturation momentum which marks the borderhne between weak and 
strong scattering, we shah now compute the structure function in the 'saturation region' at 

< (5^(x,0). This win ahow us to understand how the energy-momentum sum rule (|3^ ) IS 
satisfied and speculate about a possible partonic interpretation. 

Let us start by rewriting the scattering amplitude defined in Eq. ( |3.36| ) in terms of the 
saturation momentum introduced in the previous subsection. It reads 

T{z,x,b±) = i{l -exp [iQ'^Ql{x,b±)z'^] } (3.55) 

In this high energy, or relatively low regime, that we are interested in, there is a central 
region of the nucleus which looks 'black' to the projectile. That is, for any given z there is a 



'black disk' radius R\:){z,Q,x) such that, for all b± < i?b the eikonal phase in Eq. ( |3.55D can 
be neglected since rapidly oscillating and the scattering amplitude becomes purely imaginary 
with a magnitude equal to one, meaning that the dilaton component at radial distance z is 
completely absorbed in the shockwave. Recalling also Eq. (|3.53| ), we see that this black disk 
radius is determined by 

/(O) Q'Q'A^,0)z^' ^ ' 

For given x and Q^, a black disk exists at the center of the 'nucleus' {b± = 0) only for those 
components of the dilaton which have penetrated far enough in the radial direction, namely up 
to z > zo with zo = 1/ \/ QQs{x, 0). Physically, this means that the transverse size r_|_ ~ z of the 
partonic fluctuations of the projectile should be large enough for their color dipolar interactions 
in the target to be strong. Since on the other hand z is restricted to z < 1/Q by the dilaton 
wavefunction, it is clear that a black disk can form only when the energy is high enough for 
the condition < Q'^{x,0) to be satisfied. Once that the black disk appears at the center of 
the 'nucleus' (for a given z > zq) it rapidly occupies the whole central region of the shockwave, 
which is essentially flat in b±. When further increasing the energy, the black disk extends al 
larger values of b±, at a rate which depends upon the transverse profile of the shockwave. 
For example, for a type-II shockwave the black disk radius will be given by 

(RkAf = [Q'Q^(x, 0)z4]i/3 _ 1 . (3.57) 

This exhibits a rapid increase with the energy, proportional to (1/x)^/'^, which reflects the power 
law tail in b± of the respective energy density profile. On the contrary, a type~I shockwave with 
an exponential tail at large b± will have a black disk radius which at high energy increases with 
1/x only logarithmically, that is 

{R^Af = ln''[Q^Ql{x,0)z% (3.58) 



For high enough energy, the total cross section ( p. 42]) for a given value of z will be proportional 
to the area of the respective black disk. (Indeed, outside the black disk, at b± » Rh{z,Q,x), 
the amplitude is very small and rapidly decreasing with b±.) So, clearly, the cross-section 
associated to the type-II shockwave will rise like a power of the energy, in violation of the 
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Proissart bound, whereas that corresponding to a type-I Shockwave with exponential tail will 
saturate the Froissart bound at high enough energy. In that sense, the type-II shockwave is 
a more realistic model for a nucleus in QCD. In fact, the above mechanism for the emergence 
of the Froissart bound is exactly the same as expected in QCD: the competition between the 
power-law increase of the scattering amplitude with the energy at a given b± and the exponential 
decrease of the nuclear matter distribution at large b±. 

The previous considerations lead to the following expression for the structure function at 
Q^<Ql{x,0) : 

j\t2 ^6 rco 

i^(x,Q') = ^^y dzzKl{Qz)[R,,{z,Q,x)A]\ (3.59) 

The Bessel function will effectively cut off all contributions coming from z > l/Q, so that the 
integration in Eq. ( ^.59 ) is in practice supported only in the interval 



< ^ < A • (3.60) 



VQQ six, 0) ^ ^ Q' 

However, the dilaton wavefunction diverges strongly at small z, more precisely K^Q^;) — 
Thus for any profile f{b±A) which leads to a black disk area (R^^A)"^ larger than z"^ at 
z <C 1, therefore including the two cases in Eqs. ( p. 57 ) and ( 3.5^ ), the integration is dominated 



by its lower limit. This comes as a surprise, since one might have expected to dominant contri- 
bution to come from z ~ 1/Q, as happens in the case of the 7^~current (see Sect. ^ below), and 
also in perturbative QCD (after identifying z with the size of the dipole fluctuation of the virtual 
photon). Coming back to the case under study, we see that for z ^ zq the black disk area (RhA)'^ 
is a number of 0(1) and thus the integration is simply proportional to 1/Q^Zq = Qs{x,0)/Q^. 
So, quite remarkably, the structure function for the dilaton DIS is not sensitive to the expansion 
of the black disk with increasing energy — rather, it is controlled by the central part of the 
black disk at bj_ < 1/A — , and hence it is not affected, e.g., by the violation of the Froissart 
bound in the case of the type-II shockwave. This should be contrasted with the corresponding 
situation for the 7^-current, to be described in Sect. ^. 

As we show in Appendix |C|, deeply at saturation (Q^ <^ Ql{x,0)), we can also determine 
the precise normalization of the structure function, and we finally have 

16v2vr xAL Q 

with K a number which depends only on the profile of the shockwave under consideration. 

It is straightforward to check that in the transition region at Q ~ Qs{x), or equivalently 



X ~ Xs{Q), the above result in Eq. (3.61) is parametrically consistent with the respective extrap- 



olation of the dominant result at large-x, i.e. the 2-graviton exchange contribution in Eq. ( |3.51| ): 

XsFs{xs,Q^) ~ ~ x,f(2)(x„Q2). (3.62) 

Of course, for x ~ Xs{Q) all the terms in the twist-expansion become parametrically of the same 
order and our various approximations fail to apply there, except for parametric estimates. 
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By using these estimates, one can now verify that the energy-momentum sum rule ( 3.23| ) 
is indeed satisfied parametrically. The integral over x in the l.h.s. of Eq. ( 3.23| ) is dominated 
by a; ~ Xs{Q), since for x ^ Xg the 2-graviton exchange contribution in Eq. ( p. 51 ) is rapidly 
decreasing with x, whereas for x < Xg, we have xF{x,Q'^) ~ l/y^. Thus, by making use of 
Eqs. (|3.62|) and ( 3.54 ) one can write 



dxxF{x,Q^) ~ xiF{x„Q^ 



2 a2 



AL 



(3.63) 



Let us conclude this discussion by suggesting a possible partonic interpretation for the previ- 
ous results. The standard OPE analysis for the correlator ( |3.3[) shows that, in the weak coupling 
regime and for sufficiently high Q^, the quantity xF[x,Q'^)/Q'^ has the meaning of the gluon 
distribution in the target (since gluons are the partons which directly couple to the relevant op- 
erator J = (1/4)F^). Namely, it is proportional with the number of gluons per unit transverse 
area having longitudinal momentum fraction x and transverse momenta k± < Q. Then what 
about the previous results at strong coupling ? Since the energy-momentum sum-rule is con- 
trolled by X ~ Xs{Q), it is natural to look at Eq. ( 3.62 ), which implies XsFg{xs, Q'^)/Q'^ ~ Q'^N^. 
As we shall argue in more detail in Sect. 4^, where a similar result occurs in the context of 
the 7^-current, this estimate is indeed consistent with a partonic picture, in which the partons 
(here, gluons) are distributed along the saturation line with occupation numbers of order one 

One the other hand, there are clearly no partons at x ^ XsiQ), as obvious from the lack of 
leading-twist contributions to F{x,Q'^) in that region. 



4. DIS off the Shockwave: the 7^— current 

Let us now turn to the problem of deep inelastic scattering of an 7?.-current off the shockwave. 
As discussed in the beginning of Sect. ^ this is a conserved current in = 4 SYM which is 
carried by the fermionic and scalar fields. It is a bilinear in both of these types of fields, whose 
precise form is not necessary for our discussion (it is given, for example, in [34]), and it has mass 
dimension 3. 

4.1 General equations: from the 5D action to the polarization tensor 

We will be interested in calculating the polarization tensor 

n'^^(g) = dVd^-dVe^"'' ie{xo){[J^{b + r/2),r{b-r/2)]), (4.1) 

which is the analogous to Eq. ( p.4| ) and where the notation is the same as in the previous section. 
The bulk field dual to the 7^-current is the S0(6) gauge field A^, where we have suppressed 
the fixed color index of this vector field corresponding to the fixed index of the 7?--current. To 
compute H^'^ (see Eq. ( [4. 3] ) below), we need only the quadratic part of the respective Yang-Mills 
action in the AdSs background. It reads 

S = --^2^ I d'^xdz^gg^^g^'^FmnF^, (4.2) 

^'*At this stage, it might be useful to remember that, in the saturation region at a; < Xs(Q), the partonic inter- 
pretation of the DIS structure functions does not rely on the validity of the twist expansion (and the dominance 
of the twist-two operators) not even at weak coupling [6]. 



- 26 - 



with Fmn = dm-A-n — dnAm- According to the AdS/CFT correspondence the 2-point function of 
the current appearing in the integrand Eq. ([4.1|) is equal to 



(4.3) 



In the above, Sc\ is the action evaluated on the solution to the equations of motion 

dm(.V^r''9''''Fp<i) = 0, (4.4) 
i.e. the Maxwell equations in curved space-time, with the following boundary conditions 

limA^{z, x) = A^{x) and lirn ^^(z, x) = 0. (4-5) 

It will be convenient to adopt the gauge condition A^ = 0. Furthermore, as in the case of the 
dilaton, we shall assume that the boundary field is a plane wave of the form 

^^(x)=e-^(''"^''+''''^"U^, (4.6) 



so that the components of are pure numbers. Now let us integrate Eq. ( |4.2| ) by parts. Using 
the Maxwell equations (^^) only the boundary term survives to give 



Sci = I d^xdza„(^5"^V^i^p,^n). (4.7) 



Substitution of our shockwave metric into the above leads to 

Ar2 



'S'cl 



c 



I d*x - {-A+A'_ - A_A'^ + AA'i) _ , (4.8) 



where a prime represents differentiation with respect to z, the index i = 1,2 refers to the 
transverse components and we have dropped a term —hA-^-A'_^ inside the parenthesis which does 
not contribute since h cc for small z. 

For simplicity, from now on we will restrict to the case where the fields are independent of 
x±, which is the situation when the shockwave is homogeneous in transverse directions: dih = 0. 
Furthermore, since the shockwave metric in Eq. (|2.3D does not depend on x"*", we can restrict 
ourselves to a single Fourier mode of the gauge field as follows 

A^{z, x+,x-) = e-'^~^^A^{z, X-). (4.9) 

The Maxwell equations for n = —,z,i become 

{zd,z-^d, + iq-d-)A+ = {q~fA_, (4.10) 

(9_ - iq-h)A'^- = iq-A'_, (4.11) 

{zd,z-^d, + 2iq-d-)Ai = -{q-fhAi. (4.12) 

Notice that in the vacuum case h = 0, the x~ dependence of the field is that of a plane wave, 
that is e"*"^^^ . The first two equations are coupled. In order to solve them we differentiate 
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( [4.10 ) with respect to z and we make use of (4.11) to arrive at the differential equation which 
determines : 

{d^zd^z-^ + 2ig-9_)A'+ = -[(ffhJ^^. (4.13) 
As in the dilaton case, we separate the total field into a vacuum piece and a scattering piece 

A^{z, X-) = Af (z, X-) + Af{z, X-), (4.14) 



with A 



(0), 



Z, X 



satisfying the vacuum version of the Maxwell equations which are obtained by 



setting h = 0. We shall construct Af^{z,x ) so that its boundary condition at z = is fully 
encoded in the vacuum field, that is 



lim Af^\z,x ) = Af_,{x ) and lim ^^''^(z, a; ) = 0. 



Now we define the boundary to bulk propagator D^'^ according to 



(4.15) 



4^)(z,x-)= / dy~ D,,''{z,x-,y-)A,{y-), 



(4.16) 

which must be determined by solving the equations of motion. Then the contribution of the 
scattering piece to the classical action Eq. ( [4.8| ) is obtained from —A^_^^A'^^^—A^^^A'!f_^^ +Af^ Af^^ 
and reads 



5ci-5o = F^y dx-dy-A^{x-)A,{y-)^d,D^'''{z,x-,y-) 



(4.17) 



2=0 



with the volume factor V arising from the integration over x^ and x±. Functional differentiation 
according to Eq. ( |4.3| ) and use of the same prescription as the one used in the dilaton case (cf. the 
discussion just after Eq. ( |3.20| )) leads to^^ 

1 1 



V 167r2 z 



d,D^-{z,x-,y-) 



2 = 



Substituting into Eq. (4.1) we finally arrive at 

7 N^c 



n^-(g) 



L 167r2 



j dx-dy- e*''^^^ -y ^^d,D^"'iz,x-, 



y 



(4.18) 



(4.19) 



2=0 



Let us notice here that not all the components of the propagator D^'^ , and therefore of the 
tensor 11'^'^, are independent. In fact, only two components are independent since 11'*^ can be 
decomposed as 



ni(x,Q^) + (p'^-^< 



p 



p_q 



qnU2{x,Q^), (4.20) 



which follows from general symmetries (namely the current conservation together with the time- 
reversal symmetry Il^^'^(q) = n^''(— g)) and can be explicitly verified by calculating one by one 
all the components of the propagator D^^"^. Hence we will just need to determine H** and, for 
example, 11 , since the remaining non-vanishing components will be then determined as 



4{q- 



n- 



and 



2{q- 



n- 



(4.21) 



^^Notice that there is a factor 1/V for each functional differentiation since S/SA{x) — {l/V)5/SA{x 
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4.2 The classical solution 



The solution to the equations of motion is constructed in the same way as in the dilaton case. 
Resumming the graviton exchanges to all orders we find for the scattering part of the transverse 
field 



-2q- 



^GTiz,z;x-)T{z)AfiO,z), 
z 



(4.22) 



where the scattering amplitude T(i) has already been defined in Eq. ( |3.36| ) and with the Green's 
function Gt satisfying in momentum space 



{zd,z~^d, - K^)Gt{z, z; K^) = z6{z - z), 



(4.23) 



with the notation = —2k^q . The precise form of the Green's function is given in Appendix 
S Now we can see that the respective boundary to bulk propagator reads 



D'\z,x-,y-) 



-2q 



di 



-GT{z,z;x-)T{z)D^^{z,-y-), 



(4.24) 



where Dq is the corresponding boundary to bulk propagator in the vacuum and in momentum 
space is given by 



Di}{z,K) = KzKi{Kz). 



(4.25) 



Substituting the above into Eq. ( 4.19| ) we can immediately perform the integrations over the 
longitudinal coordinates which will select the component of both the Green's function Gt 
and the propagator Dq. We then have (recall that = —2q^q^ > 0) 



167r2 xAL 



^-dzGT{z,z;Q'^ 
z z 



2 = 



T{z)D'o\z,Q). 



(4.26) 



Finally by making use of 



1 



d,GT{z,z;Q^) 



z=0 



QzKiiQz) 



(4.27) 



and Eq. (4.25), and changing the name of the integration variable from z to z, we arrive at 



167r2 xAL 



dzzKi{Qz)T{z). 



(4.28) 



We recognize here a structure similar to that in Eq. (|3.4lD : once again, the Bessel function 
Kf^Qz) plays the role of the wavefunction squared for the partonic fluctuation of the 7?.-current 
(namely, a fluctuation with transverse size Ax± ~ z). Interestingly, it is the same Bessel function 
which determines the wavefunction of the quark-antiquark fluctuation of a space-like photon in 
lowest-order perturbative QCD. The same remark applies to the structure of IT , that we now 
compute. 

Namely, the scattering part of the component of the gauge field is found as 



Af\x-,z) 



-2q- 



di 



^GL(z,i;x-)r(i)4°)(0,i), 



(4.29) 
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with the Green's function Gl satisfying in momentum space 

id,zd,z'^ - K^)Gl{z, z; K^) = z5{z - i), 



(4.30) 



and again the precise form of the Green's function is given in Appendix |^. Now we need to be 
careful when we express the vacuum field yli°^(0 , z) in terms of the boundary values. We have 

A!!^\x-,z) = I dy-d,Do''{z,x-,y-)A^{y-), (4.31) 

where the upper index /x can be either + or — . Both components of the boundary to bulk 
propagator dzDQ^{z,x~ ^y~) satisfy the vacuum version of Eq. (4.13) and in momentum space 
they are given by 

2{q-f 



1 



dzD^^{z,K) = -K^zKq{Kz)c^' with c+ = 1 and c 



(4.32) 



The coefficients have been determined by plugging the above equations into the left hand side 
of Eq. ([4.101) and then evaluating both sides of that equation on the boundary. Since we would 
like to calculate 11 we shall need here only the Dq component, but it is already obvious 
from the above equation why the ratio 11"' /IT is given by the second equation in Eq. ( 4.2l| ). 
Now we are ready to give the boundary to bulk propagator which reads 



di 



[z,x ,y ) = -2q 
We proceed as before and we plug into Eq. (l4.19[) to obtain 



-Gl{z,z;x )T{z)diDQ {z,-y ). 



n 



167r2 xAL 



1 

-^-GUz,£;Q') T{z)d,Di}{z,Q). 

Z Z 2=0 



Using —{1/z)Gi^{z,z;Q'^)\z=o = zKq{Qz) and then letting z ^ z we finally obtain 



n- 



167r2 xAL 



dzzKl{Qz)T(z). 



(4.33) 



(4.34) 



(4.35) 



4.3 Structure functions and their partonic interpretation 

We are finally in a position to calculate the Lorentz invariant structure functions Fi and F2. In 
general they are given by 

Fi = ^lmTli, F2 = -^lmU2. (4.36) 
ZTT zvr 

Focusing on the ( ) and (ii) components in ( [4.20| ) we can express Hi and 112, and therefore the 

structure functions, in terms of the components 11 and 11** that we have already calculated. It 
is also customary and convenient to introduce the transverse and longitudinal structure functions 
i^T and Fl. Let us summarize here our results for all these various structure functions: 

Ft = 2xFi = - Imn- = ^ ^ J dz z K^Qz) T{z), (4.37) 
Fl = ^Im^n- = -JL-^ dzzKl{Qz)T{z), (4.38) 



vr {q-y 167r3 AL 

Ft + Fl. (4.39) 
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These results are in agreement with those in Ref. [40]; there is in fact a mismatch by an overah 
factor AL, but this is presumably related to the use of different conventions in the normalization 
of the current-current correlator. 

We have already seen that the scattering amplitude T that emerged during the calculation 
is the same as the one appeared in the dilaton case. Thus the saturation scale will also be 
the same in the case we are currently considering. This is natural since the saturation scale 
is a property of the nuclear target and hence it should be the same whether we probe it with 
the gluon current or with the 7^-current. Recalling that = ttA^L/2x (in a homogeneous 
situation) we can write the amplitude T as 

r(z,Q) =i[l-exp(iQ2Q2^4^j _ 

Since it not possible to obtain an analytic expression for the structure functions for arbitrary 
value of X and Q^, we shall separately consider the two limiting cases. 

(i) For Q » Qs{^) (or, equivalently, x » Xs{Q), with Xs{Q) defined in Eq. ( |3.54] )) we need 
to expand the exponential in (4.40) to second order since the first order term will not contribute 



to the imaginary part of the tensor H^'^ . The upper limit in the z-integration is l/\^QQs 3> 1/Q, 
so it can be set equal to oo. We find 

''- = ^-'$, ^- = ^-^ F, = ^-^. (4.41) 

As in the dilaton case, these contributions are of higher-twist order (in the context of the TZ- 
current, a leading-twist contribution would be independent of at large Q^, up to logarithms), 
which reflects the absence of point-like constituents in the nucleus. 

(ii) For Q <C Qs{x) (or, equivalently, x <C Xs{Q)), we set the exponential in ([4.40| ) equal 
to zero and we restrict the integration to z > l/\/QQs- In the longitudinal sector, and since Kq 
diverges only logarithmically at small z, the integration is dominated by the region z ~ 1/Q, so 
that the lower limit can be safely set equal to zero. This is not the case in the transverse sector 
where the integration is sensitive to small values. Since Ki{Qz) ~ 1/Qz for small z, there is a 
logarithmically enhanced contribution when compared to i*L. Therefore the transverse sector 
dominates in F2, albeit this dominance is only logarithmic. We get 

Remarkably, the structure functions at low have the same parametric form as for a proton in 
QCD at weak coupling, except for the replacement of the baryonic color factor Nc corresponding 
to a proton by a factor A'^^ corresponding to the density of degrees of freedom in the shockwave. 
Even though the above expressions in Eqs. ( 4.411) and ( 4.42 ) are strictly valid for Q ^ Q 



s 



and Q <^ Qs respectively, they become of the same order when extrapolated to Q ~ and 
therefore they can be parametrically trusted even for Q around Qs- This also means that for 
Q ^ Qs the transverse and longitudinal structure functions are of the same order. 

So far, our analysis in the 7^-current case has been done for a homogeneous, in the transverse 
space, Shockwave. However, given the close analogy of the structure functions in Eqs. ( [4. 371) and 
( [4.38 ) to the dilaton structure function in Eq. ( |3.41| ) and the common form of the scattering 
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amplitude in the two cases, it is not unreasonable to extend the discussion of Sect. 3.4 to the 
problem of an inhomogeneous shockwave probed by an 7^-current. In such a situation, and 
according to Eq. ( 3.59| ), we would write the transverse and longitudinal structure functions at 



saturation (Q^ <^ Q'^(x,0)) as 

i^T,L(x,Q') = ^ J dzzKlo{Qz)[RUz,Q,x)A]'', (4.43) 



with i?b and zq as in Sect. 3.4. But, in contrast to the dilaton case, the z-integration is now 
dominated by 2 ~ 1/Q in the longitudinal sector, while it is only logarithmically sensitive to 
the lower limit in the transverse sector (like in the homogeneous case), a behavior similar to the 
one in QCD. For the exponential profile leading to the black disk radius in Eq. ( |3.5^ ), we find 

where we see again the extra logarithmic enhancement in the transverse structure function. 
In what follows, we shall argue that the results in Eq. ( 4.42| ) have a natural physical in- 



terpretation in terms of parton saturation in the boundary gauge theory. Our discussion will 
closely follow that of the plasma structure functions at strong coupling, as originally presented 
in Ref. [34]. To that purpose, two more ingredients turn out to be extremely useful: 

(a) The associated energy-momentum sum rules, which can be derived via the same strategy 
as in the dilaton case (recall the discussion of Eq. ( |3.50| ) ) . Namely, one first computes the leading- 
twist, or single scattering, contributions to the Lorentz invariants Iii(x, Q"^) and I\.2{x,Q'^)-, which 
are real quantities as anticipated: 

n^P^^^ a„d n?> = il^^, (4.45) 

Then one uses analyticity and contour integration in the complex u = \/x plane to finally deduce 

dxxFi(x,Q2) = _^A2 and j ^ dx F^ix , Q^) = K\ (4.46) 

which should more properly be understood as limiting values when Q"^ ^ co (since, as in 
the dilaton case, we have neglected the higher-twist contributions to Il{v,Q^) near v = 0). 
Examining Eqs. (4.41) and ( [4.42| ), we can verify that these sum rules are parametrically satisfied 



by our previous estimates for the structure functions, and that the integrals are dominated by 
X ~ Xs-, as expected. In fact, by using the exact, integral, expressions for 11** and 11 derived 
in the previous subsection, one can explicitly evaluate the above integrals of Fi and F2 and thus 
check not only their high-Q^ limit shown in Eq. ( [4.46 ), but also the subleading corrections at 
large Q^. This will be detailed in Appendix 



Once again, the sum-rules ( 4.46 ) extract the contribution of the energy-momentum tensor 
to the operator product expansion for the current-current correlator. In view of this and of the 
mass dimensions of the above integrals, it is natural to interpret them as the energy per unit 
length in the nucleus rest frame. 
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(b) The expression of the ghion distribution at saturation in perturbative QCD, which 
reads, parametricahy [6] 



- 02 In '^^^^^ 



(4.47) 



where we now assume A = g^N^ <^ 1 and Nc ^ 1. (The last assumption is not essential, but 
merely convenient for the sake of comparison with the AdS/CFT results.) In writing Eq. ( 4.47| ), 
we have also divided by the nuclear transverse area ttR'^, so that the quantity shown there is 
the number of gluons per unit transverse area having transverse momenta k± < Q < Qs{x). 
Then, the estimate in the r.h.s. of this equation can be understood as follows: the gluons with 
longitudinal momentum fraction x and low transverse momenta k± < Qs{x) are saturated, in 
the sense that there is a number of gluons of 0(1/A) of each color per unit cell in the phase- 
space, and this number does not grow with the energy anymore (in contrast to what happens 
at larger momenta kx_ > Qs{x))- Specifically, Eq. ( [4.47| ) implies the following estimate for the 
gluon occupation number at saturation in QCD at weak coupling and large Nc : 



n(x, kj 



1 



dxG{x,Q^) 



dQ2 



1 , Q|(x) 
A fc2 



(4.48) 



which shows only a weak, logarithmic, dependence upon the energy, via the corresponding 
dependence of the saturation momentum: Ql{x) ~ with cj ~ 0(A). 



We shall now argue that the gluon distribution in Eq. ( 4.47 ) is the natural quantity to 
compare with the structure function F2 in Eq. ( 4.42| ) and that this comparison suggests a 
partonic interpretation for the latter which is further supported by the sum rules O). 

As already mentioned, the 7^-current directly couples to fermionic and scalar fields of = 
4 SYM, so the respective structure function F2 is most naturally related (at least at weak 
coupling, where the standard OPE arguments apply) to the parton distributions for these fields. 
However, the latter lie in the adjoint representation of the color group; hence, for the purposes 
of the Nc power counting, it is more appropriate to compare F2 va. M = A SYM to the gluon 
distribution in QCD, and not to the respective structure function, which rather describes the 
distribution of quarks in the fundamental representation. Such a comparison, together with the 
formal similarity between Eq. (4.47) and the expression ( [4.42| ) for F2 (they both show the same 
parametric dependencies upon N"^ and Q"^)-, immediately suggests that partons are saturated at 
strong coupling too, but with occupation numbers of order one, rather than 1/A. More precisely, 



for k±_ <C Qs{x), Eq. (14.421) implies 



n(x 



ki.) 



AL dF2(x,Q2) 



dQ2 



In 



Q=k_L 



ki 



(4.49) 



which suggests that, at strong coupling, the mechanism for parton saturation has the same 
effect as at weak coupling [6] (although its dynamical nature might be very different) — it 
limits the rate for parton emission with increasing rapidity y = ln(l/x) : dn/dY ~ const, for 
k±_ < QsiY). At least at weak coupling, this limitation is associated with the formation of a 
'color glass condensate' [6] — a strong color mean field collectively created by all the small-x 
partons, which blocks the radiation of new partons once the density is high enough. It is likely 
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that a similar picture holds at strong coupling as well. But unlike at weak coupling, where this 
saturation requires parametrically large occupation numbers n ~ 1/A ^ 1, to compensate for 
the weakness of the coupling, at strong coupling it occurs already for occupation numbers of 
since the mutual repulsion becomes strong as soon as the partons are allowed to interact 
with each other. 

In writing Eq. ( [4.49 ) we have extracted the parton occupation number at saturation out of 



the structure function F2 in the same way as we would do at weak coupling. To further justify 
this identification, let us now show that it is also consistent with the sum rules ( 4.46| ), which 



in fact explain the additional factor AL that we have inserted in the definition of n{x^k±_) in 



Eq. (4.49). We have already mentioned below Eq. (4.46) that the following integral 



C dxF2{x,Q^) ~ XsF2{xs,Q^) ~ NlK^, (4.50) 

JO 

is proportional to the energy density per unit longitudinal distance in the target rest frame. 
Using this, we would like to construct the energy per unit transverse area in the shockwave^^. 
This is obtained by multiplying the quantity above by the longitudinal extent L (to get the total 
energy), then dividing it by the transverse area {i.e. multiplying it by a factor A^), and finally 
multiplying it by the Lorentz factor 7 (to make a boost to the target infinite momentum frame). 
One thus gets 

d^su 



-iLK^ (xF2(x,g2))^^^^ ^ k+{KLF2{x,Q^))^^^^ , (4.51) 



where we have also identified the longitudinal momentum = xp~^ = X7A of the parton (scalar 
or adjoint fermion) which absorbs the virtual photon in this frame. The last estimate above 
suggests that, at least for x ~ where most of the energy is located, the quantity ALF2{x, Q"^) 
has the meaning of the total number of partons (with energy fraction x) per unit transverse area 
as probed on the resolution scale Q^: 

ALF2{x,Q^) ~ / ^'k^^Ij^- (4-52) 



This is indeed the same as Eq. (4.49) since, by definition, n{x,k_L) = (l/N^) [dn/d^k±d^b±]. 
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A. From black holes to Shockwaves 



Let us modify the temperature dependent piece of the AdSs black-brane metric by an arbitrary 
function depending on the coordinate. That is, using Fefferman-Graham coordinates, we 
consider the metric in 

ds2 = J_/_(l^^dt2 + (l + a) rdxi + (dx3)21 +dz2\ with a = c{x^)^. (A.l) 
z"^ [ 1 + a J 4zq 

Here c{x^) is a dimensionless function, which reduces to unity for the black-brane solution, and 
zq = l/vrT with T the Hawking temperature. In general Eq. ( [A.l| ) is not a solution to the 
Einstein equations. Now we would like to make a boost and choose a frame which moves 
with velocity —v with respect to the frame x"* and along their common axis in the 3-direction. 
In the new frame we would also like to switch to light-cone coordinates and it is straightforward 
to check that the combined transformation of coordinates reads 

t = ^ x+ + 2(i±!^lf-, (A.2) 

x3=^i x+-2^x-. (A.3) 

Clearly the transverse and radial coordinates remain unaffected. The metric tensor transforms 
according to 

S-rrin = -qTH Qpq, (A.4) 

SO that the components which are transformed read 



7^(1 + , ^ _ 2 7^(1 + ^) ^0 



[gtt + 933)= I ' (-^-6) 



g+- =9-+ =2 - 933) = -;2(r^- (A-7) 

When the boost is ultra relativistic, that is for j ^ I 4^ v ^ 1, we expect the metric to simplify. 
Indeed we have 

, , [c(i-)]2z8 

n ' 1 2^'^C{X-)Z^ 1 J ^ . 1 + 16^g o^ 

•Z -Zn 1-1- ^'•^ 1^ 2: 1 -I- '^v^ r 

where we have substituted a with its explicit form and we have defined c(x~) = c(x'^). Notice 
that in this ultra relativistic limit Eq. (|A.3D simplifies to x'^ ~ —^/2'yx~. Furthermore, close to 
the boundary, i.e. when z <C 2^01 one can approximate gu ~ l/z^ in Eq. ( [A.l| ) and also neglect 



the last factor in g and g^ in Eq. ( |A.8| ). However, notice that g does not vanish since it 

is enhanced by the square of the large boost factor. Putting everything together, and dropping 
the primes, we can write the metric as 



ds^ ~ 



^^'^^^ {dx-f - 2dx+dx- + dxl + d^2 
4 



(A.9) 
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This is formally the same as the metric corresponding to a homogeneous in the transverse space 
Shockwave, cf. Eq. ( |2.6| ), which we recall is an exact solution to the Einstein equations. Here, 
however, this metric emerges only as a near-boundary approximation valid for z <^ 1/T. On 
the other hand, the more general boosted-plasma metric in Eq. (A. 8) keeps trace of the black 
hole horizon at z = zq, and it is valid only for z < zq. This makes it natural to supply the 
Shockwave in Eq. (|2.6| ) with a radial cutoff at finite z 1/T or in general at 1/A, where A is an 
infrared cutoff in the dual gauge theory. 

Let us finally observe that when the function c{x'^) has support only in an interval of finite 
length L, in the boosted system it will be transformed into a delta function. For instance, in 
the case of a slice of a plasma (in the dual gauge theory) with c{x'^) = 1 for < L/2 and zero 
otherwise, in the boosted frame and in light cone coordinates we find 

c{x-)c^^^5{x-). (A.IO) 
B. The bulk-to bulk propagator 

In this Appendix we outline the construction of the Green's functions for the equations of motion. 
For definiteness, let us consider the dilaton case. Taking the Fourier transform of ( 3.25| ) w.r.t. x~ 
and x± we have 

.2 3 



di - -d,- K^]G(z,z;K^) = z''6(z-z), (B.l) 
z ) 

where we have defined the variable = /c^ — 2q~k~^. The solution to the eigenvalue problem 
for the operator 5^ — (3/z) dz is given by 

(^d'^z - ^8^^"^ h{ujz) = -oj^z^ 32icoz), (B.2) 

and the eigenfunctions satisfy the completeness relation 

dtouj 2;^J2(a;z) i^J2(a;i) = z'^6{z — z). (B.3) 

Therefore we can write the Green's function as 

j-OO 1 

Giz,z;K^) = - dujio^ —.z'^J2{ioz)z'^32{uJz). (B.4) 

Now we would like to integrate over oj and we need to distinguish two possibilities according to 
the sign of K^. In the spacelike case > we find 

i; >{)) = -z'^f[K2{Kz)l2{Kz)Q{z - z) + K2{Kz)l2{Kz)Q{z - z)]. (B.5) 

In the timelike case < we need a prescription to avoid hitting the pole which is located on 
the positive w-axis. Letting — ie with e > 0, we find (in agreement with the analytic 

continuation of ([B.5|)) 



2 '2 

G{z, z; <0) = -'-^-^[R'^\\K\z)h{\K\z)e{z - z) + R^^\\K\z)h{\K\z)@{z - i)].(B.6) 
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Following the same procedure we can construct the Green's functions that we encountered for 
the 7^-current in Sect. ^. The transverse and longitudinal Green's functions in equations ( 4.23| ) 
and ( [4.30 ) are given by 



Gt{z,z;K^ > 0) 
Gl{z,z;K^ >0) 



-zz[Ki{Kz)Ii{Kz)Q{z -z)+ Ki{Kz)Ii{Kz)Q{z - £)], (B.7) 
-zz[Ko{Kz)lo{Kz)e{z -z)+ Ko{Kz)lo{Kz)e{z - z)], (B.8) 



respectively. The corresponding functions for < are obtained via analytic continuation. 
Now the Green's functions in coordinate space will be given by 



G{z,z;x -y ,x±-y±) 



(27r)3 



.ik+(x -y )+ik±ix^-y±)G{z,z;K^). (B.9) 



and the integration over k~^ can be easily performed using the representation in Eq. (B.4) with 
our prescription — ie. The pole is at {uJ^ + kY)/2q~ — ie and is enclosed by the 

integration contour only for > V - Thus the fc^-integration leads to the factor 



d)fc+ Q-ik+{x -y ) 



27r cj2 + A;2 - 2q-k+ - ie 



- y ) 
2q- 



exp 



i{x -y + 
2^^ 



, (B.IO) 



which in turn yields the expression in Eq. (3.26). 

The bulk-to-bulk propagator has been already constructed in the past (see for example [16, 
56,57]), but mostly using Euclidean signature, or using standard Minkowski coordinates instead 
of light cone coordinates. Therefore, let us now show how the methods used in this paper lead to 
these known results. We will aim to write the propagator as a function of the chordal distance 
introduced earlier in Sect. |2|. This property arises from the fact that the differential operator 
can be rewritten as the Laplacian of the hyperbolic space M^+i with coordinates (z, xi, . . . , Xd). 
We shall be interested in the case d = 4. Defining the 4~dimensional vector r = x — y and 
performing a Wick's rotation in Eq. (B^), the Euclidean Green's function in coordinate space 
is obtained as (below, fc^ = Yli=i is strictly positive) 



G{z,z;r) 



2 '2 

-Z Z 



dU 

(2^ 



ik-r 



duj U! 



1 



u;2 + A;2 



J2(a;z)J2(u;i) 



(B.ll) 



First we perform the integration over k. In d-dimensions we have 



dU 



Ak-r 



{2ttY Lo^ + k^ (27r)'^/2 \r 



Using the above for d = 4 and the identity Ki(a;r) 
coordinate space propagator as 



(- 

V r 



Kd/2-iiuJr). (B.12) 
■(l/a;)(9rKo(wr), we can write the 



G{z,z;r) 



z^z d 







4^2 J, Qj. 

where we have found convenient to define 



du;a;J2(w^;) J2('^-2)Ko(u;r) 



z z d (r2 — riY 



4^2 J, Qj. f-^f2{r2 + ri)2 



(B.13) 



ri 



Yr2 + (z^^i)2 and r2 = \/ r"^ + {z + i)^. 



(B.14) 
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Now we introduce a chordal coordinate ^, which is closely related to the chordal distance q 
defined in Sect. |2|, according to 

2zz 1 

^= ^2 + ^2 + ^2 = ^ ^^^^ 

It is straightforward to make a change of variable from r to ^ and performing the required 
algebraic manipulations we finally arrive at 

where 2-^1 is the hypergeometric function. As announced, this expression depends on the vari- 
ables z, z and r only through a single chordal coordinate. Up to a minus sign, which is due to 
different conventions, this is equal to the bulk-to-bulk propagator given in [57]^''. 

Furthermore, let us choose the 'transverse' space to be 2-dimensional. By making use of 
Eq. (B.12) for d = 2, we can write the propagator as 



2-2 roo 2-2 (j.^ _ 

G2{z,z;r) = -—- / d^ u;J2(u;z)J2M)KoK) = ^ ^ (B.17) 

We can easily express the above in terms of the chordal distance q to find 

zz 

G2{z,z;q)=--—^2Fi{3,5/2,5;-l/q). (B.18) 

Multiplying the above by 8tt'^E6{x^)/N^, as dictated from Eqs. ( |2.18| ) and ( |2.4D , we are lead to 
Eq. ( 2.19| ). Finally, yet another definition is to use coshu = 1/^, which means that we can write 



' —97/ 

zz e 



G2{z,z;u) = -—-—, (B.19) 
47r smhn 



which, up to an overall factor independent of n, is the form used in [16] . 



C. Prefactor of the structure function 

In this short Appendix we give an alternative calculation of the structure function at saturation 
in the limit <C Q'^{x, 0), in the case of probing with a dilaton, and in order to determine the 
correct prefactor. We can write the most dominant contribution to the structure function as 

Ar'2 r)6 A fR roo 

i^(x,Q') = ^^y^ d^b^y^ dzzKliQz){l-cos[Q^Q^,ix,b^)z'']}, (C.l) 



where the radius R is defined by Qj{x,R) = Q . We have already seen in Sect. |3^ that the 



dominant contribution in this dilaton case comes from the region z <^ l/Q, thus it is natural to 

^^One needs to set A = d = 4 in Eq. (6.12) of that reference in order to obtain the 5-dimensional massless case. 
^*In Eq. (25) in that reference, where the sign in the exponent should be negative. 
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approximate the Bessel function by its lowest order term, that is K2{Qz) ~ 2/{Qz)'^. Then the 
2;-integration gives \/27r Qs{x,b±)/Q^ and the structure function becomes 



16\/2^ xAL Q TT Jo V /(O) 



Notice that we have replaced the radius R in the upper limit of the impact parameter integration 
by oo. This happens because the presence of the profile function in the numerator of the 
integrand effectively restricts the integration at values bj_ < 1/A, while at the same time R ^ 
1/A. A posteriori this also justifies the small-z approximation since for all values of b± that 



contribute, one has <C Q'l{x,b±). Now, k in Eq. (C^) is a pure number which depends on 
the profile of the shockwave. For instance, k = 8 for the exponential type-I shockwave of the 
form f{b±A) = (l/27r) exp(— 6_lA), while k = 2 for the type-II one. 



D. Explicit check of sum rules 

Here we will explicitly check the validity of the sum rule for the structure function Fi as stated in 
the first expression in Eq. (4.46). One can follow a similar procedure for other sum rules. Using 
Eqs. (4.37) and (4.40) which provide us with Fi and the scattering amplitude T respectively we 
have 



327r3 xAL 



dz z KliQz) [1 - cos (Q2Q^ j _ 



(D.l) 



We wish to integrate the above from a; = to x = 1 which is the physical region. Making a 
change of variables from z to C = Qz and from x to t = Q1{x)/Q'^ , we can write 



dx xFi{x, 



647r 



dCCK?(C) 



oo 

dt — 



cos tC, 



t2 



(D.2) 



with ti = Q2(i)/q2 ^ T^^t^LjlQ'^. It is straightforward to see that the sum rule will be exactly 
recovered only when we set ti — > as the lower limit in the t-integration. For small but non- 
vanishing ti, we can separate an integration from to cxd which gives the sum rule, minus an 
integration from to ti <C 1 for which we can do a Taylor expansion. We find 



dx xFi(x, 



NlA^ mlk^ QliX) 



+ 



(D.3) 
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